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PREFACE 


x  Four  short  presentations  were  made  to  the  Sub-Committee  on  Aeroelastieity,  during 
the  50th  Meeting  of  the  Structures  and  Materials  Panel,  on  the  results  of  a  cooperative 
programme  on  active  flutter  suppression  on  a  dynamic  mode)  of  the  YF-17  aircraft. 

During  this  programme,  British  Aerospace,  MBB,  Northrop,  Air  Force  Flight  Dynamics 
Laboratory,  NASA  and  ONERA  cooperated  in  deriving  control  laws  for  active  flutter 
suppression  for  one  explosive  wing-store  flutter  case  of  the  model.  Phase  control  laws  were 
all  tested  and  compared  during  wind  tunnel  tests  performed  in  the  Langley  16  ft  wind  tunnel 
in  October  1979.  Results  were  quite  promising  and  open  the  way  for  future  cooperation  on 
full-scale  aircraft.  -  ;  -  ■  s  V  , 

The  Sub-Committee  on  Aeroelastieity  unanimously  agreed  that  the  four  presentations 
should  be  published  as  an  AGARD  Report  since  they  were  thought  to  be  an  important 
contribution  to  the  development  of  active  control  technology. 


G.COUPRY 

Chairman,  Sub-Committee 
on  Aeroelastieity 
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WIND  TUNNEL  TEST  OF  A  FIGHTER  AIRCRAFT  WING/STORE 
FLUTTER  SUPPRESSION  SYSTEM  -  AN  INTERNATIONAL  EFFORT 

C.  Hwang,  E.  H.  Johnson,  G.  R.  Mills 
Northrop  Corporation 
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Hawthorne,  California  90250 

T.  E.  Noll 

Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  Air  Force  8ase,  Ohio  45433 

M.  G.  Fanner 

NASA  Langley  Research  Center 
Hampton,  Virginia  23665 


SUMMARY 

A  30%  scale,  half  span  model  of  a  lightweight  fighter  aircraft  with  an  active  wing/store  flutter  sup¬ 
pression  system  was  tested  in  the  NASA  Langley  Research  Center  Sixteen  Foot  Transonic  Dynamics  Tunnel  in 
September-October  1979.  The  test  featured  a  store  configuration  that  was  intentionally  designed  to  exhibit 
a  violent  flutter  condition.  In  addition  to  Northrop  organized  control  laws,  three  European  countries  also 
contributed  control  laws  to  stabilize  this  condition.  After  the  control  laws  were  mechanized  by  Northrop, 
they  were  tested  at  the  Langley  facility.  The  model  was  tested  up  to  170%  of  the  open  loop  flutter  dynamic 
pressure  in  a  number  of  cases,  with  the  indication  that  a  substantially  greater  improvement  was  achievable. 
This  paper  discusses  some  special  features  of  the  test  model,  the  design  and  implementation  of  the  control 
laws,  the  test  monitoring  techniques  and  detailed  test  results. 

INTRODUCTION 

Wing/store  flutter  problems  continue  to  degrade  the  performance  of  fighter-attack  aircraft  with  air- 
to-ground  combat  roles.  Flutter  speed  placards  caused  by  the  carriage  of  many  combinations  of  external 
store  configurations  have  reduced  aircraft  speed  capability  and  as  a  result,  decreased  its  survivability. 
Active  flutter  suppression  has  shown  promise  of  preventing  store  flutter  and  has  the  potential  for  provid¬ 
ing  significant  improvements  in  aircraft  operational  and  mission  effectiveness.  Since  flutter  suppression 
systems  can  use  electro-hydraulic  feedback  networks  to  counteract  the  flutter  conditions,  the  eventual 
integration  of  the  suppression  system  into  an  advanced  aircraft  control  system  holds  promise  of  extending 
the  aircraft  performance  with  a  minimal  increase  in  control  hardware. 

In  order  to  develop  a  working  suppression  system  dealing  with  wing/store  flutter  and  to  resolve  tech¬ 
nical  problems  in  actual  systems  implementation,  an  active  flutter  suppression  system  was  designed  and 
incorporated  into  a  30%  scale  model  of  the  YF-17  aircraft.  The  model  has  a  leading  edge  and  a  trailing 
edge  control  surface.  A  store  configuration  was  intentionally  designed  that  featured  a  violent  flutter 
condition.  In  the  series  of  wind  tunnel  tests  performed  in  1977,  a  feedback  control  system,  consisting  of 
a  number  of  filters  and  a  simple  gain,  phase  shift  circuit  was  able  to  suppress  violent  flutter  conditions. 
For  instance,  at  M  =  0.60,  the  projected  improvement  in  the  flutter  dynamic  pressure  was  29%.  The  analyt¬ 
ical  and  test  results  related  to  the  test  Sequence  are  documented  in  an  AGARD  report  (Reference  1). 

The  experience  gained  in  the  above-mentioned  wind  tunnel  tests  pointed  the  way  to  further  improve¬ 
ments  that  could  be  made  in  the  flutter  suppression  system  and  the  test  apparatus:  1)  Augmenting  and 
upgrading  the  mechanical  and  electrical  components,  2)  designing  mote  complex  control  laws  with  the 
goal  of  a  demonstrated  improvement  of  70%  in  the  flutter  dynamic  pressure,  and  3)  using  the  Hewlett- 
Packard  5451C  Fourier  Analyzer  to  monitor  the  control  law  performance  and  to  point  out  ways  they  can 
be  improved. 

In  order  to  concentrate  the  new  test  effort  on  flutter  control  law  development  and  suppression  system 
performance,  a  single  test  configuration  was  selected  which  featured  the  above-mentioned  severe  flutter 
condition.  In  addition  to  Northrop-developed  control  laws,  three  control  laws  were  supplied  by  European 
countries  at  the  invitation  of  the  U.S.  Air  Force.  The  international  contributors  were  British  Aerospace, 
MB8,  and  0NERA.  These  control  laws  were  interpreted,  mechanized  and  integrated  into  the  flutter  sup¬ 
pression  electronics  console  at  Northrop's  Hawthorne  facility  prior  to  test  entry. 

The  wind  tunnel  tests  were  performed  in  September-October  1979  at  the  NASA  Langley  Center  Sixteen- 
Foot  Transonic  Dynamics  Tunnel.  In  addition  to  Northrop  and  AFFDL  personnel,  all  European  participants 
were  on-site  to  witness  the  tests  and  to  suggest  modifications  to  their  control  laws,  as  the  need  arose. 
This  paper  describes  the  major  features  of  the  improved  flutter  suppression  model,  the  test  monitoring 
techniques,  and  the  analytical  and  experimental  results  of  the  newly-developed  control  laws.  The  control 
laws  developed  by  the  European  participants  are  presented  separately  by  their  originators. 

THE  REFURBISHED  MODEL  AND  NEW  INSTRUMENTATION  FEATURES 

The  flutter  suppression  model  is  a  30%  scale  half  span  model  mounted  on  a  set  of  roll  bars  and  stabi¬ 
lized  with  forward  and  aft  tension  cables.  Leading  edge  and  trailing  edge  control  surfaces  are  operated  by 
miniaturized  hydraulic  actuators.  A  detailed  arrangement  of  the  model  is  shown  in  Figure  1  with  additional 
information  on  the  model  given  in  References  1,  2. 

The  basic  configuration  adopted  in  the  present  test  sequence  is  Configuration  B  of  the  previous  test 
sequence.  It  features  an  empty  tip  launcher  rail  and  an  AIM-7S  missile  attached  to  a  pylon  near  the  wing 
tip  (W.S.  60.75).  As  part  of  the  refurbishing  work,  the  leading  edge  and  trailing  edge  surface  hinge  mech¬ 
anisms  were  completely  redesiqned  to  improve  their  ruggedness.  The  original  angular  potentiometers  attached 
to  the  hinges  were  replaced  with  rotary  variable  differential  transformers  (RVTD's).  The  new  design  with 


RVDT's  achieved  a  more  reliable  surface  position  sensing.  It  also  featured  new  surface-to-actuator-shaft 
attachments  for  improved  reliability  and  interchangeability. 


Four  (4)  Sundstrand  accelerometers  were  installed  in  the  wing  sections  to  serve  as  the  sensors 
for  flutter  suppression.  The  sensors,  identified  as  a,  through  a.,  were  installed  in  the  following 
locations:  1  ' 


Accelerometer  No. 

1  (ax) 

w.s. 

51.45 

F.S. 

145.18 

(25%c) 

Accelerometer  No. 

2  (a2) 

w.s. 

51.45 

F.S. 

158.00 

( 76%c ) 

Accelerometer  No. 

3  (a3) 

w.s. 

60.75 

F.S. 

148.56 

( 25%c ) 

Accelerometer  No. 

4  (a4) 

w.s. 

60.75 

F.S. 

157.47 

( 701c ) 

The  approximate  locations  of  these  sensors  are  indicated  in  Figure  1.  Additional  model  instrumen¬ 
tation  included  wing  and  horizontal  tail  moment  strain  gages,  the  control  surface  servo  inputs,  and  the 
angular  position  potentiometers  mentioned  above.  For  the  new  test  sequence,  additional  wiring  and  tub¬ 
ing  were  installed  to  accomnodate  the  wing  and  tail  moment  gage  bridges  and  the  instrumentation  of  the 
AIM-7S  store  functioning  as  the  flutter  stopper. 

The  control  console  panel  was  modified  to  accommodate  additional  required  controls  and  indicators. 
The  existing  patch  boards  were  reorganized  and  the  control  console  was  rewired  to  accommodate  new  control 
laws  developed  at  Northrop  and  by  European  country  participants.  Additionally,  in  order  to  improve  the 
actuator  response,  a  new  hydraulic  power  unit  was  acquired  to  supply  fluid  at  the  rate  of  1.6  gallons  per 
minute  at  a  pressure  of  2000  psi. 

Remote  control  of  the  hydraulic  power  unit  was  added  to  the  console  to  provide  for  start/stop  and 
bypass  control  of  this  unit  from  the  wind  tunnel  control  room. 

THE  FLUTTER  DETECTOR  AND  THE  FLUTTER  STOPPER 

Two  pieces  of  hardware  were  added  to  the  flutter  suppression  system  in  order  to  protect  the  model  in 
case  of  an  unexpected  flutter  occurrence.  The  first  of  these  was  a  flutter  detector  which  was  designed  to 
indicate  when  a  flutter  condition  existed,  while  the  second  device  was  a  "flutter  stopper"  which  featured 
a  mass  inside  the  AIM-7S  store  that  could  be  fired  ahead  rapidly  to  interrupt  the  flutter  condition  once 
it  was  identified  by  the  flutter  detector.  Separate  bays  in  the  control  console  monitored  and  controlled 
the  performance  of  the  flutter  detector  and  the  flutter  stopper. 

The  flutter  detector  has  a  number  of  electronic  circuits  which  sense  when  the  amplitude  of  the  in¬ 
put  signals  from  the  model  are  equal  to  or  above  a  preset  level.  A  digital  counter  counts  the  number  of 
pulses  detected.  A  triggering  signal  is  generated  when  a  preset  number  of  pulses  are  accumulated  by  the 
detector.  To  prevent  false  triggering  due  to  spurious  noise  signals,  the  pulse  counter  is  reset  to  zero 
if  a  certain  time  lapse  occurs  between  consecutive  pulses.  In  wind  tunnel  tests,  when  the  trigger  condi¬ 
tions  are  met,  the  triggering  signal  may  be  applied  to  activate  the  flutter  stopper  and  to  operate  the 
tunnel  flow  bypass  valves.  Figure  2  is  a  block  diagram  showing  the  basic  operation  of  the  flutter  detec¬ 
tor.  Figure  3  is  the  oscillograph  record  of  a  number  of  the  responses  of  the  YF-17  model  instruments 
acquired  in  the  wind  tunnel  when  a  moderate  flutter  condition  was  encountered.  The  time  traces  read 
from  left  to  right,  with  all  transducers  and  scales  identified  at  the  left  side  of  the  figure.  The  top 
trace  of  the  figure  shows  the  trigger  signal,  which  appeared  only  after  the  wing  bending  and  torsion 
bridge  gages  experienced  a  number  of  high  stress  cycles.  Prior  to  that  time,  the  stress  cycles  with  a 
moderate  amplitude  were  ignored  by  the  flutter  detector. 

A  scale  model  of  the  AIM-7S  Sparrow  store  was  designed  to  function  as  the  flutter  stopper.  The  new 
store  featured  a  spring-loaded  movable  mass  (weight  *  3  pounds).  In  the  normal  latched  position,  the 
store  model  has  the  same  c.g.  location  and  mass  inertia  in  pitch  and  yaw  as  the  store  used  in  the  previous 
tests  (Reference  1).  In  case  of  severe  flutter,  the  mass  can  be  released  by  a  solenoid  mechanism,  which 
allows  the  mass  to  move  forward  inside  the  store  housing  for  a  maximum  distance  of  10  inches  in  approxi¬ 
mately  1,  second.  The  store  with  the  changed  c.g.  location  was  designed  to  disrupt  and  stop  the  flutter 
condition  of  Configuration  B.  When  the  movable  mass  is  in  the  forward  position,  it  can  be  reset  remotely 
to  any  desired  latched  position  along  the  10-inch  travel  distance.  The  reset  is  executed  through  the 
operation  of  a  gear  motor  and  ball-nut  assembly.  A  linear  potentiometer  with  a  contact  point  on  the 
movable  mass  indicates  the  instantaneous  mass  location  inside  the  flutter  stopper.  A  sectional  drawing 
showing  the  flutter  stopper/guide  rail  assenbly  is  given  in  Figure  4.  Figure  5  shows  the  details  of  the 
assembly  inside  the  store. 

Flutter  analyses  of  the  model  were  performed  when  the  mass  was  in  the  normal  and  deployed  positions. 
Figure  6  shows  the  results  of  these  analvses.  It  Indicates  that  the  flutter  stopper  can  increase  the 
model  flutter  speed  quite  siqni ficantly.  Two  different  curves  are  shown  for  the  deployed  result.  In  one, 
the  analysis  was  performed  without  any  modification  to  the  analytical  data,  while  in  the  other,  the  bend¬ 
ing  and  torsion  frequencies  for  the  system  with  the  mass  deployed  were  acquired  in  a  ground  vibration  test 
and  input  to  override  the  analytical ly-obtained  frequencies.  It  is  seen  that  the  increase  in  flutter 
speed  with  the  measured  frequencies  is  less  beneficial  than  with  the  purely  analytical  result.  Neverthe¬ 
less,  both  results  provide  for  a  satisfactory  pressure  margin  in  flutter  dynamics.  The  effectiveness  of 
the  flutter  stopper  was  later  confirmed  in  the  wind  tunnel  entry. 

NORTHROP  fONTROl  LAW  DEFINITION  AND  ANALYSIS 

In  order  to  establish  accurate  open  loop  aircraft  transfer  functions  as  a  basis  for  new  control  law 
development  and  to  make  use  of  the  previous  wind  tunnel  test  data,  a  now  flutter  analysis  was  performed  on 
the  YF-17  flutter  suppression  mode).  Simultaneously,  a  ground  vibration  test  was  performed  at  t tie 
Hawthorne  facility.  During  the  test,  the  tunnel  supporting  system  used  for  the  YF-17  model  was  duplicated 
to  the  extent  possible.  Table  1  presents  a  comparison  among  the  natural  vibration  frequencies  measured 


previously  in  the  ground  vibration  tests  at  Wright-Patterson  Air  Force  Base  and  at  the  Hawthorne  facility, 
and  those  computed  in  the  present  analysis.  This  analysis  differs  slightly  from  the  final  analysis  pre¬ 
sented  in  Appendix  A  of  Reference  1  in  that  the  stiffnesses  used  in  that  analysis  had  been  arbitrarily 
reduced  slightly  in  order  to  match  the  test  results.  It  is  seen  that  the  present  analysis  shows  a 
reasonable  correlation  with  the  test  results. 

TABLE  1.  NATURAL  VIBRATION  FREQUENCIES  (IN  Hz)  OBTAINED  IN  THE  GROUND  VIBRA¬ 
TION  TESTS  COMPARED  WITH  ANALYTICAL  RESULTS,  CONFIGURATION  B. 


MODE 

AFFOL  GVT 

HAWTHORNE  GVT 

TUNNEL  GVT 

ANALYSIS 

1 

4.80 

4.44 

4.40 

4.63 

2 

7.40 

7.10 

7.40 

7.19 

3 

12.80 

12.31 

12.20 

13.40 

4 

14.20 

14.23 

14.00 

14.88 

5 

19.00 

18.66 

19.40 

19.14 

A  study  was  performed  to  determine  the  level  of  refinement  necessary  in  the  wing  paneling  for 

unsteady  aerodynamic  forces  in  order  to  obtain  a  converged  solution.  Special  attention  was  given  to  the 

proper  paneling  of  the  two  control  surfaces  used  in  flutter  suppression.  Based  on  this  study,  a  panel 
scheme  using  195  boxes  on  the  wing  surface  was  selected  for  the  analysis.  The  new  analysis  for  Configura¬ 
tion  B  resulted  in  a  matched  point  dynamic  pressure  of  84  psf  at  M  =  0.80  (see  Figure  6).  The  flutter 
frequency  at  the  match  point  is  6.0  Hz. 

Once  confidence  in  the  open  loop  analysis  had  been  attained,  the  construction  of  new  control  laws 
could  proceed.  Figure  7  is  a  representative  block  diagram  of  the  new  control  laws.  In  the  figure,  6-.  and 
t>Q2  are  two  model  response  outputs  (e.g.,  the  angular  and  vertical  accelerations  derived  from  the 
outboard  sensor  pair)  which  are  filtered,  blended  and  then  fed  back  to  one  of  the  control  surfaces.  G,  and 

G2  are  aircraft  transfer  functions  that  related  model  response  to  the  command  to  the  control  surface 

actuators.  The  n  terms  represent  response  to  noise  sources,  primarily  the  tunnel  turbulence.  The  F^  terms 
are  pseudo-integrators  which  convert  the  accelerometer  outputs  to  responses  which  approximate  velocity  and 
displacement  signals.  The  "pseudo"  designation  is  used  because  the  integrators  are  of  the  form  l/(s+a) 
where  "a"  is  a  constant  used  to  prevent  steady  state  drift.  The  R,  components  are  gain  settings  which  con¬ 
trol  the  relative  amplitudes  of  the  velocity  and  displacement  signals,  with  the  R5  term  controlling  the 
overall  gain.  The  H  block  is  a  series  of  filters  which  act  together  as  a  bandpass  to  limit  the  amount  of 
signal  which  is  fed  back  at  frequencies  removed  from  the  flutter  frequency. 

A  synthesis  procedure  (Reference  3)  was  applied  to  determine  the  constants  Rj  through  R4.  Alto¬ 
gether,  four  control  laws  were  organized  by  Northrop  and  tested  in  the  wind  tunnel.  They  are  identified  as 
control  laws  Nl,  N3,  N3T  and  N3P.  The  control  surfaces  and  the  accelerometers  used  by  each  of  these  con¬ 
trol  laws  and  the  control  laws  contributed  by  the  European  countries  are  identified  in  Table  2. 


TABLE  2.  CONTROL  LAWS  FOR  THE  FLUTTER  SUPPRESSION  MODEL 


CONTROL 

LAW 

RESPONSIBLE 

ORGANIZATION 

CONTROL 

SURFACE 

ACCELEROMETERS 

LEADING 

EDGE 

TRAILING 

EDGE 

D 

a2 

a3 

a4 

Nl 

No rthrop 

l 

H 

M 

N3 

II 

9  jj 

N3T 

" 

1 

X 

H 

■ 

N3P 

•1 

■ 

HI 

BAE 

British  Aero¬ 
space 

R 

D 

hBB 

Messerschmi tt- 
Bol kow-Binhm 

D 

S 

ON 

0NERA 

I 

■ 

Q 

The  block  diagrams  of  the  four  Northrop  control  laws  are  shown  in  Figures  8  through  11. 

Pretest  analyses  were  performed  on  all  Northrop  control  laws,  with  some  additional  analyses  made  at 
Hawthorne  while  the  test  was  in  progress  at  the  NASA  Langley  TDT  facility.  Figures  12  through  15  present 
analytical  data  for  the  four  Northrop  control  laws.  Figure  12  is  a  plot  of  the  stability  behavior  of  the 
open  loop  system  and  of  Control  Law  Nl  at  M  -  0.80  using  dynamic  pressure  Q  as  a  parameter.  The  data  were 
acquired  by  the  characteristic  diagram  method.  Figures  13  a,  b  are  the  open  loop  Nyquist  plots  of  Control 
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Law  N3  at  two  dynamic  pressures.  Based  on  Nyquist  data  of  this  type,  the  system  gain  and  phase  margins 
were  established  as  a  function  of  the  dynamic  pressure.  Figures  14  a,  b  are  similar  open  loop  Nyquist 
plots  for  Control  Law  N3T.  The  data  Indicated  the  N3T  Control  Law  was  substantially  less  effective  than 
Law  N3.  This  point  was  confirmed  In  the  wind  tunnel  tests  and  Is  discussed  later  In  the  paper.  Control 
Law  N3P  used  both  the  leading  and  trailing  edge  control  surfaces  and  is  simply  a  superposition  of  the  N3 
and  N3T  control  laws.  The  stability  plot  of  this  law,  based  on  the  characteristic  diagrams  method,  is 
shown  in  Figure  15. 

CONTRIBUTED  CONTROL  LAWS 

As  mentioned  in  the  previous  section,  three  (3)  control  laws  were  contributed  by  three  separate 
European  country  participants:  The  British  Aerospace  Corporation,  MBB  and  ONERA.  The  control  laws  are 
identified  by  the  letters  BAE,  MBB  and  ON,  respectively  (see  Table  2).  The  block  diagrams  representing 
the  three  control  laws  are  given  in  Figure  16  through  18.  As  can  be  readily  observed,  even  though  they 
were  designed  to  suppress  a  single  flutter  condition,  the  control  laws  feature  a  great  diversity  in  their 
basic  organizations  and  in  the  transducers  and  control  sui  face(s)  they  used. 

TEST  MONITORING  TECHNIQUES 

An  important  aspect  of  the  project  was  the  development  of  procedures  that  monitored  and  directed  the 
wind  tunnel  test  program.  Prior  to  the  test  entry,  monitoring  programs  were  developed  on  the  Hewlett- 
Packard  5451C  Fourier  Analyzer.  Actual  servo  hardware  and  an  analog  simulation  of  the  YF-17  model  were 
used  to  try  out  and  debug  the  monitoring  programs.  In  this  section,  some  key  monitoring  programs  used 
in  the  tunnel  entry  are  described. 

Nyquist  Plot 

A  key  monitoring  program  was  developed  to  extract  the  open  loop  transfer  function  of  the  aircraft 
flutter  suppression  system  when  the  control  loop  was  closed.  The  following  sketch  defines  the  terminology 
that  is  used  in  the  discussion. 


where 


■  actuator  and  servo  transfer  function 
G  »  open  loop  transfer  function  of  the  aircraft 
H  *  transfer  function  of  the  fepdhack  loop 


A  further  term  is  Gc =  the  closed  loop  transfer  function: 


G 


gag 

c  '  1  +  GaGH 


(1) 


The  Nyquist  criterion  required  a  plot  of  G^GH.  The  assumption  for  developing  the  test  monitoring 
technique  was  that  it  should  be  possible  to  make  a  relatively  accurate  measurement  of  Gc  but  that  a  direct 
measurement  of  G  would  be  difficult  since  it  would  require  the  measurement  of  a  mode  that  was  suppressed 
and  because  noise  would  contaminate  a  measurement  made  inside  the  loop.  The  alternative  was  to  extract  G 
from  the  measured  values  of  Gc,  G ft  and  H.  (The  latter  two  transfer  functions  could  be  measured  very  well 
with  the  loop  open.)  The  open  loop  transfer  function  was  then  obtained  from  Equation  1: 


G  = 


GA  '  GcH> 


(2) 


The  Nyquist  diagram  could  then  be  calculated  directly. 


The  above  example  illustrates  the  case  involving  a  single  feedback  signal.  For  a  multiple  sensor 
system,  reference  is  made  to  a  typical  flutter  suppression  control  law  as  shown  in  Figure  7.  In  the 
figure,  the  F^'s  and  H  are  compensation  filters,  while  the  R^’s  are  potentiometer  settings.  In  this  case, 
the  measurements  are  made  of  Gjc  and  G2C,  which  are  6qj/<5c  and  002/<5c>  respectively.  The  open  loop  trans¬ 
fer  functions  Gj  and  G2  are  then: 


r  _ _ _ ulc  _ _ _ _ 

1  '  r^TH  [G1c  (R-|  F1  +  R2  Fz)  ♦  C,2c  "(R3  F3  Vr4  F4)i 

_ p2,c _ 

G2  1  •  R5H  [G,c  (Rf  f,  +  r?‘f2)  +  g2c  (r3  f3  +  r4  f4TT 


(3) 

(4) 


Once  the  open  loop  transfer  functions  have  been  computed,  it  is  possible  to  construct  the  open  loop 
Nyquist  plot.  The  formula  for  the  Nyquist  function  is: 


N  =  R5H  [ G-j  (R1  F1  +  Rp  Fr>)  +  {^3  ^3  +  ^4  F4)  1 


(5) 
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The  application  of  the  Nyquist  plot  technique  in  actual  wind  tunnel  test  will  be  illustrated 
later  in  the  section  dealing  with  Experimental  Results. 

The  measurement  of  the  closed  loop  transfer  functions  required  that  some  type  of  excitation  be 
applied  at  the  6C  point  of  Figure  7.  During  the  simulation  work  carried  out  at  Hawthorne,  as  well  as 
during  the  later  wind  tunnel  test,  the  excitation  utilized  was  a  random  transient  that  was  generated 
by  the  Fourier  Analyzer.  This  was  a  random  noise  signal  that  was  set  to  zero  in  the  last  25  percent 
of  the  sample  period.  This  made  the  response  to  the  excitation  die  out  in  the  sample  period  ( i . e . ,  it 
was  periodic  in  the  "window")  and  fast  Fourier  transform  techniques  were  valid.  Four  thousand  samples 
were  taken  for  every  sample  period  at  a  rate  of  200  samples/sec.  for  a  total  of  20  seconds  per  period. 

The  maximum  frequency  measured  with  this  setup  was  100  Hz.  The  number  of  averages  used  varied  depend¬ 
ing  on  the  number  required  to  give  an  adequate  signal  to  noise  response.  In  later  wind  tunnel  tests, 
a  minimum  of  six  (6)  ensemble  averages  were  used. 

Smoothing  of  Transfer  Functions 

Another  test  monitoring  program  dealt  with  the  smoothing  of  frequency  domain  signals  such  as  aircraft 
transfer  functions.  Because  of  the  substantial  level  of  turbulence  in  the  tunnel,  the  response  signal  to 
noise  ratios  are  such  that  key  modal  information  may  be  masked  by  the  superimposed  noise.  In  order  to 
reduce  the  noise  effect,  smoothing  techniques  were  applied  to  the  measured  or  processed  transfer  functions 
Methods  of  smoothing  entail  performing  a  Fourier  transform  on  the  transfer  function  to  obtain  the  impluse 
response  of  the  system.  It  is  assumed  that  all  the  response  after  a  certain  time  in  this  record  is  uncor¬ 
related  with  the  input  excitation.  A  simple  means  of  removing  this  repsonse  is  to  apply  a  rectangular 
window  to  the  impulse  response  that  zeroes  out  all  response  after  a  certain  time. 

When  this  truncated  response  is  transformed  back  to  the  frequency  domain,  the  resulting  transfer 
function  is  considerably  cleaner,  but  retains  all  the  key  modal  information. 

A  criticism  of  this  "rectangular"  window  is  that  it  creates  a  high  frequency  signal  that  is  superim¬ 
posed  on  the  transfer  function.  This  behavior  can  be  eliminated  if  a  different  type  of  window,  the 
exponential  window,  is  applied  to  the  impulse  response. 

A  drawback  of  an  exponential  window  is  that  it  affects  the  estimated  damping.  Reference  4  shows 
that  if  this  window  is  applied  to  a  single  degree  of  freedom  transfer  function,  the  damping  is  modified 
according  to  the  formula: 


c'/c 


o 


c/c0 


(6) 


Where  c'/c0  is  the  damping  obtained  from  the  smoothed  transfer  function,  c/c0  is  the  actual  damping, 
n1  is  the  decay  constant  in  the  exponential  window  and  uin  is  the  natural  frequency  of  interest.  If  the 
window  is  applied  n  times,  the  correction  term  is  also  increased  by  a  factor  of  ' n ’ .  Although  this  does 
not  appear  to  be  aunajor  correction,  it  was  decided  to  avoid  this  complication  and  use  the  rectangular 
window. 

Peak-Hold  Damping  Trend 

The  peak-hold  spectrum  method  provides  Fourier  transformed  measurements  of  the  model  responses,  fil¬ 
tered  through  250  narrow-band  circuits.  The  peak  response  within  each  interval  is  registered  on  a  screen 
until  the  peak  amplitudes  are  stabilized.  The  damping  of  a  resonant  mode  is  proportional  to  the  inverse 
of  the  peak  amplitude  at  that  mode  frequency.  The  damping  trend  data  can  be  extrapolated  as  the  dynamic 
pressure  Q  is  increased  to  determine  the  projected  flutter  dynamic  pressure  for  either  a  controlled  or 
uncontrolled  system.  Because  of  the  simplicity  of  the  approach  in  that  knowledge  of  the  excitation  is 
not  required,  the  method  was  used  extensively  in  the  tunnel  test  to  determine  the  control  system  effec¬ 
tiveness  and,  with  the  assistance  of  other  monitoring  programs,  to  revise  the  tunnel  test  strategy. 

Other  monitoring  programs  developed  for  the  wind  tunnel  tests  included  the  determination  of  new 
feedback  constants  to  cause  a  moderate  phase  shift  of  the  feedback  signal,  and  an  optimization  program 
to  improve  the  control  signal  while  the  test  was  in  progress.  Additionally,  open  and  closed  loop  trans¬ 
fer  function  acquisition,  as  well  as  other  processing  techniques,  were  applied  to  the  real  time  data 
recorded  on  a  14  channel  FM  tape  system.  A  substantial  portion  of  the  results  presented  in  the  following 
section  was  generated  from  the  FM  tape  data. 

EXPERIMENTAL  RESULTS 

Peak-Hold  Damping  Data 

After  setting  up  the  flutter  suppression  model  in  the  NASA  Langley  TDT  facility  and  performing  a  sim¬ 
ple  ground  vibration  test,  actual  wind-on  tests  were  initiated  in  early  October  and  lasted  15  working 
days.  Altogether,  27  runs  were  made  consisting  of  511  test  points.  The  key  results  for  the  uncontrolled 
flutter  model  and  for  the  model  with  various  control  laws  are  surrmarized  in  Figures  19  through  26  in  the 
form  of  peak-hold  damping  trend  plots.  The  plots  were  based  on  the  signals  from  the  wing  torsion  moment 
strain  gage. 

Referring  to  Figure  19,  the  damping  trend  indicated  a  projected  open  loop  flutter  dynamic  pressure  of 
75  psf.  This  result  was  confirmed  by  two  other  methods.  One  method  was  based  on  the  Nyquist  criteria 
using  open  loop  data  extracted  from  the  closed  loop  transfer  functions  as  explained  in  the  previous  sec¬ 
tion  on  Test  Monitoring  Technqiues.  The  other  was  the  Royal  Aircraft  Establishment  (RAE)  method  based  on 
the  extrapolation  of  the  two  natural  modes  in  flutter.  Experimental  results  throughout  the  test  period 
showed  a  consistent  flutter  dynamic  pressure  of  75  psf.  Post  processing  with  the  HP  modal  analysis  indi¬ 
cated  a  flutter  dynamic  pressure  of  74  psf.  This  was  in  contrast  to  the  predicted  analytical  result  of 
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84  psf,  and  the  experimental  results  of  83.5  psf  acquired  in  1977.  It  is  difficult  to  determine  the  pre¬ 
cise  cause  of  the  deviation  in  the  flutter  dynamic  pressures.  The  most  probable  cause  was  that  the  model 
dynamic  behavior  was  sensitive  to  the  constraints  on  the  half  fuselage.  The  constraints  were  caused  by 
the  supporting  system  which  was  somewhat  redundant.  Any  slight  misalignment  could  influence  the  natural 
frequency  modes  in  a  subtle  manner,  which  in  turn  could  affect  the  flutter  dynamic  pressure.  Another  fac¬ 
tor  which  could  affect  the  analytical  results  were  the  deviations  in  physical  data  of  the  AIM-7S  Store. 
After  test,  slight  deviations  were  discovered  of  the  new  store  mass,  c.g.  and  pitch  inertia  properties. 
Using  the  corrected  data,  the  flutter  dynamic  pressure  was  dropped  from  84  psf  to  80  psf.  During  the  test 
entry,  because  of  test  time  limitations,  and  in  view  of  the  consistent  experimental  results,  it  was  decided 
to  bypass  this  problem  and  to  use  a  dynamic  pressure  ratio,  (Q/Qf),  for  evaluation  purpose  when  analytical 
and  experimental  data  were  compared. 

Figures  20  through  26  show  that  all  the  control  laws  were  effective  in  suppressing  the  flutter  condi¬ 
tion.  On  the  figures,  the  symbols  indicate  actual  data  points  while  the  straight  lines  are  projections 
obtained  from  a  linear  regression  of  the  data  points. 

RMS  Response  Data 

The  FM  tape  data  was  processed  to  calculate  power  spectra  information  for  the  majority  of  the  test 
points  where  peak  hold  information  was  taken.  The  purpose  of  this  analysis  was  to  provide  a  more  quantita¬ 
tive  comparison  of  the  various  control  laws.  It  also  provides  insight  into  the  model  behavior  and  into 
subjective  criteria  that  were  used  to  guide  the  conduct  of  the  test. 

For  the  leading  edge  control  laws.  Table  3  presents  RMS  values  for  the  leading  edge  control  surface 
position  and  the  wing  root  bending  and  torsion  strain  gage  data.  The  results  are  presented  in  two  fre¬ 
quency  ranges,  allowing  for  the  distinction  between  the  lower  frequency  response  associated  with  the  open 
loop  flutter  mode  and  those  responses  associated  with  higher  frequencies.  The  first  row  in  Table  3  gives 
the  responses  at  the  highest  open  loop  dynamic  pressure  at  which  the  model  was  tested.  A  careful  examina¬ 
tion  of  the  data  in  this  table  reveals  some  of  the  subtleties  of  the  various  control  laws.  For  instance, 
it  is  seen  that  the  N1  and  the  BAE  control  laws  did  a  better  job  of  suppressing  the  flutter  mode  than  the 
N3  law,  but  that  the  N3  law  had  considerably  less  response  in  the  10-20  Hz  range  than  the  aforementioned 
laws.  Figures  27-29  give  the  power  spectral  responses  of  the  leading  edge  surface  position  for  the  three 
control  laws  and  make  this  point  even  more  clearly.  A  final  point  from  Table  3  is  that  the  open  loop 
response  of  Test  Point  64  is  quite  high  and,  in  fact,  indicates  that  this  was  not  a  prudent  condition  at 
which  to  be  testing. 

Figures  30  and  31  compare  the  PSD's  of  the  wing  torsion  gage  for  the  uncontrolled  case  of  Test  Point 
64  with  the  controlled  case  at  a  dynamic  pressure  that  is  75%  higher.  Note  that  there  is  20  db  magnitude 
difference  in  the  scales  of  the  two  figures.  The  uncontrolled  PSD  is  dominated  by  the  response  at  6  Hz, 
while  for  the  controlled  case,  the  6  Hz  mode  is  suppressed  so  well  that  higher  frequency  modes  dominate 
the  response. 

TABLE  3.  RMS  RESPONSE  DATA  OF  THE  YF-17  MODEL  FOR  VARIOUS  LEADING  EDGE  CONTROL  LAWS. 


CONTROL  LAW 

Q 

(PSF) 

■  |  v 

LEADING  EDGE  POSI¬ 
TION  (DEGREES) 

WING  BENDING  MOMENT 
(IN-LBS) 

WING  TORSION  MOMENT 
(IN-LBS) 

0-10  Hz 

0-20  Hz 

0-10  Hz 

0-20  Hz 

0-10  Hz 

0-20  Hz 

None 

73 

64 

- 

- 

41.4 

49.8 

42.7 

44.1 

N1 

98 

97 

.52 

.86 

29.9 

51.0 

8.3 

22.1 

N3 

101 

124 

.63 

.69 

37.6 

50.2 

11.7 

17.7 

BAE 

105 

201 

.61 

.96 

29.9 

41.6 

8.3 

22.6 

ON 

100 

471 

.96 

.96 

62.0 

70.1 

16.5 

22.6 

N1 

128 

104 

.70 

1.10 

41.7 

61.0 

9.0 

28.6 

N3 

128 

127 

.87 

.97 

49.5 

66.2 

12.5 

22. 1 

BAE 

129 

211 

.71 

1.13 

37.3 

51.5 

9.6 

27.7 

Table  4  shows  RMS  data  for  the  trailing  edge  control  laws,  and  the  N3P  law,  which  used  both  control 
surfaces.  The  trailing  edge  laws  are  seen  to  have  required  less  control  activity,  but  they  were  not  able 
to  demonstrate  as  much  of  an  improvement  in  dynamic  pressure  as  the  leading  edge  laws.  Transfer  function 
data  tend  to  indicate  that  the  trailing  edge  laws  were  characterized  by  low  phase  margins,  making  them 
sensitive  to  small  changes  in  the  feedback  law  or  the  tunnel  condition. 

The  N3P  control  law  performed  well,  particularly  when  it  is  considered  that  no  attempt  was  made  to 
improve  upon  the  law  during  the  tunnel  test.  The  leading  edge  RMS  response  for  this  law  is  given  in  Table  4  where 
it  is  seen  to  be  considerably  less  than  for  any  of  the  single  surface  laws  at  an  equivalent  test  condition.  Figure 
32  shows  the  PSD  of  this  response  which,  when  compared  with  Figures  27  and  28,  shows  that  the  peak  ampli¬ 
tudes  are  much  less  for  the  N3P  law,  although  there  is  response  over  a  broader  frequency  range. 

Unfortunately,  trailing  edge  data  was  not  recorded  onto  FM  tape  for  the  N3P  law,  but  the  impression 
formed  during  the  test  was  that  the  trailing  edge  amplitude  was  less  during  the  test  of  the  N3P  law  than 
it  was  for  any  of  the  trailing  edge  laws  performing  alone. 
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TA8LE  4.  RMS  RESPONSE  DATA  OF  THE  FLUTTER  SUPPRESSION  MODEL  FOR  VARIOUS  TRAILING  EDGE  CONTROL  LAWS  AND  THE  N3P  LAW. 


CONTROL  LAW 

TEST 

POINT 

CONTROL  SURFACE 
INPUT  (DEGREES) 

WING  BENDING  MOMENT 
(IN-LBS) 

WING  TORSION  MOMENT 
(IN-LBS) 

0-10  Hz 

0-20  Hz 

0-10  Hz 

0-20  Hz 

0-10  Hz 

0-20  Hz 

N3T 

96 

249 

.32 

.34 

28.2 

39.0 

26.2 

29.8 

MBB 

99 

356 

.24 

.25 

45.0 

53.0 

28.0 

30.2 

N3P 

128 

N3P 

.52 

.66 

42.0 

59.0 

18.0 

26.9 

A  final  comment  is  that  the  PSD  responses  measured  from  the  FM  tape  could  have  been  measured  equally 
well  during  the  wind  tunnel  test.  This  would  have  complemented  the  peak  hold  measurements  that  were  used 
to  decide  whether  the  test  could  proceed  to  higher  dynamic  pressures.  The  HP  Fourier  Analyzer  has  consid¬ 
erable  flexibility  in  terms  of  available  displays  and  calculations  it  can  perform.  Data  gathered  during 
the  test  could  be  stored  on  disk  and  recalled  for  display  as  needed,  thereby  centralizing  and  streamlining 
the  real  time  data  gathering  process. 

Closed  Loop  Transfer  Function  Data  of  the  YF-17  Model 

During  the  wind  tunnel  tests,  model  transfer  functions  were  acquired  through  control  surface  excita¬ 
tion.  The  excitation  of  either  the  leading  edge  or  the  trailing  edge  control  surface  was  initiated  by  the 
random  DAC  signal  from  the  Hewlett-Packard  Fourier  Analyzer.  The  responses  were  measured  by  the  accelerom¬ 
eters  and  the  wing  root  section  strain  gage  bridges.  The  DAC  and  the  response  signals  were  recorded  on  the 
FM  tape  for  later  processing.  Simultaneously,  selected  transfer  functions  were  created  using  six  ensemble 
averages  of  20  seconds  each  during  the  tests  and  were  stored  in  the  Hewlett-Packard  System  Disc.  The 
transfer  functions  generated  onsite  were  used  to  monitor  the  model  behavior,  and  for  a  closed  loop  system, 
to  determine  the  effectiveness  of  the  control  law.  To  differentiate  from  the  on  site  generated  transfer 
functions,  the  transfer  functions  generated  at  Hawthorne  based  on  the  FM  data  (and  using  the  same  number 
of  ensemble  averages)  are  called  the  post-test  transfer  functions.  The  transfer  functions  presented  in 
this  section  were  acquired  with  Control  Law  N3P  in  operation,  i.e.,  the  control  law  using  both  the  leading 
and  trailing  edge  control  surfaces.  The  responses  used  in  generating  the  transfer  functions  were  the 
acceleration  differentials  (83  -  aj),  U4  -  83),  as  well  as  the  wing  root  section  bending  moment  and  tor¬ 
sional  moment  gages.  These  data  corresponding  to  Q =  95  psf  are  presented  in  Figures  33  through  38. 

Data  in  Figures  33  through  36  were  generated  onsite,  while  Figures  37  and  38  are  data  from  post-test 
processing. 

Referring  to  the  figures,  the  major  response  peaks  below  20  Hz  were  higher  using  the  trailing  edge 
surface  than  with  the  leading  edge  surface  (by  a  factor  of  two  or  more).  The  difference  is  most  prominent 
for  transfer  functions  involving  (a3  -  ai).  (Compare  Figures  33,  35.)  Also,  the  phasing  signal  is 
substantially  noisier  when  leading  edge  excitations  are  used. 

The  effect  of  increasing  dynamic  pressure  Q  for  the  same  transfer  function,  is  that  the  peak  ampli¬ 
tudes  vary  approximately  with  the  dynamic  pressure.  (Data  not  presented  in  this  paper.)  For  a  set  of 
transfer  functions,  the  changes  of  the  phase  signal  corresponding  to  moderate  Q  changes  (from  72  to  95  psf) 
were  barely  noticeable. 

An  examination  of  the  transfer  functions  for  the  wing  bending  moment  and  torsional  moment  gages 
reveals  generally  similar  patterns  for  both,  some  mixing  of  the  actual  bending  and  torsional  sig¬ 
nals  by  the  strain  gage  bridges.  As  to  be  expected,  the  modal  peaks  in  the  12-15  Hz  range,  corresponding 
to  the  pylon-yaw/wing-bending  mode,  and  the  tail  surface  bending  mode,  were  more  prominent  for  the  transfer 
functions  of  the  bending  moment  gage.  (See  Figure  37.) 

Extraction  of  Open  Loop  Response 

The  need  to  smooth  the  measured  model  transfer  functions  prior  to  other  processing  was  discussed  pre¬ 
viously.  The  noisiness  of  the  transfer  functions  may  be  evidenced  in  the  data  presented  in  Figures  33 
through  38.  This  subsection  illustrates  the  extraction  of  the  open  loop  transfer  functions  from  the  closed 
loop  transfer  functions  after  they  are  smoothed  through  the  use  of  the  rectangular  filters. 

Closed  loop  transfer  functions  were  acquired  for  Test  Point  504  at  M  =  0.80,  Q  =  95  psf,  with  Control 
Law  N3  in  force.  After  smoothing,  the  results  are  given  in  Figures  39,  40.  In  order  to  extract  the  open 
loop  transfer  functions  for  Nyquist  criteria  analysis,  it  is  necessary  to  establish  the  feedback  functions 
F-j H ( i  =  1,4)  as  shown  in  the  block  diagram  of  Figure  7.  The  four  functions  are  presented  in  Figure  41, 
while  the  Rj  parameters  for  this  case  are,  respectively,  0.365,  0.116,  0.351,  0.837  and  0.750.  Figures 
42  and  43  show  two  open  loop  transfer  functions  that  result  when  all  these  data  are  entered  into  Equations 
3  and  4. 

With  the  open  loop  transfer  functions,  it  is  possible  to  construct  the  open  loop  Nyquist  plot  from 
which,  for  the  case  of  a  single  control  surface,  the  system  gain  and  phase  margins  can  be  readily  deter¬ 
mined.  The  Nyquist  function  is  generated  according  to  Equation  (5).  For  the  case  described  above  (Figures 
39-43),  the  resulting  Nyquist  plot  is  shown  in  Figure  44.  It  is  seen  that  the  plot  has  the  desirable 
encirclement  around  the  -1  point,  indicating  Control  Law  N3  at  this  test  condition  has  large  gain  and  phase 
margins. 

Durinq  the  wind  tunnel  test,  a  number  of  Nyquist  plots  were  acquired  in  the  above  described  manner. 
They  were  valuable  in  determining  the  margins  of  the  control  system,  as  well  as  in  defining  ways  of 
improving  control  laws. 
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Determination  of  System  Frequencies  and  Dampings 

The  previous  sub-sections  describe  how  the  measured  transfer  function  data  are  prepared  for  identifi¬ 
cation  of  the  system  frequency  and  damping  parameters.  The  HP  Fourier  Analyzer  has  a  modal  analysis  pack¬ 
age  (Reference  5)  that  is  well  suited  for  the  identification  task.  One  drawback  of  that  system  is  that  it 
cannot  identify  modes  that  are  negatively  damped.  Hewlett  Packard  hasprovided  Northrop  with  a  modified 
version  of  the  system  that  does  have  a  capability  for  identifying  unstable  modes. 

Another  complicating  factor  in  the  identification  process  is  that,  near  the  flutter  point,  there  are 
two  lightly  damped  modes  that  are  very  close  in  frequency.  This  can  confuse  the  identification  system, 
with  significant  user  interaction  required  to  obtain  satisfactory  results.  With  these  qualifications. 
Figure  45  shows  open  loop  frequency  and  damping  estimates  from  the  data  gathered  during  Run  #27  of  the 
tunnel  test.  Only  the  two  lowest  frequency  modes,  which  represent  most  of  the  total  flutter  mechanism, 
are  shown.  The  dynamic  pressure  range  is  limited  by  the  range  over  which  it  was  prudent  to  excite  the 
model.  The  flutter  Q  is  74  psf  for  this  case  with  a  rapid  crossing  of  the  stability  boundary. 

Figure  46  compares  the  open  loop  result  with  those  obtained  when  Control  Law  N3  (see  Table  2)  was 
used.  It  is  seen  that  over  the  range  measured,  the  system  remains  very  stable  and  in  fact  it  is  impos¬ 
sible  to  project  a  closed  loop  flutter  speed.  The  comparison  of  the  open  and  closed  loop  frequencies 
shows  that  the  modes  remain  well  separated  when  the  control  law  is  in  force. 

Figure  47  shows  a  similar  comparison  for  Control  Law  N3T,  a  trailing  edge  law.  No  closed  loop  data 
are  available  for  this  control  law  above  80  psf,  but  the  implication  is  that  while  this  system  is  still 
stable,  it  does  not'perform  as  well  as  the  N3  law. 

Finally,  Figure  48  shows  the  results  when  the  N3  and  N3T  laws  are  engaged  simultaneously.  This  is 
designated  the  N3P  law,  and  it  is  seen  that  this  law  gives  the  best  stability  characteristics.  In  fact 
the  modal  damping  is  so  high  that  identification  is  considerably  less  precise  for  this  case. 

COMPARISON  OF  ANALYTICAL  AND  EMPIRICAL  TRANSFER  FUNCTIONS 

The  development  of  the  control  laws  relied  upon  the  analytically  derived  transfer  functions  that 
relate  the  accelerometer  responses  to  the  control  surface  inputs.  It  is  of  interest,  therefore,  to  see 
how  closely  transfer  functions  measured  in  the  tunnel  compare  with  the  analytical  ones.  A  confounding 
factor  is  that  the  measured  open  loop  flutter  dynamic  pressure  was  approximately  12%  lower  than  the  pre¬ 
dicted  flutter  dynamic  pressure,  a  fact  mentioned  previously  in  the  paper.  The  results  given  below  are, 
therefore,  presented  as  a  function  of  Q/Qf,  where  Q  is  the  dynamic  pressure  at  which  the  transfer  function 
is  presented  and  Qf  is  the  open  loop  flutter  dynamic  pressure  (predicted  or  measured,  depending  on  the 
context) . 

Figures  49,  50  present  the  comparison  between  the  measured  and  analytical  open  loop  transfer  func¬ 
tions  of  the  (aj  -  ai)  response  due  to  the  leading  edge  input.  The  comparison  was  made  at  a  (Q/Qf)  ratio 
of  approximately  1.30.  The  experimental  open  loop  data  were  extracted  from  the  closed  loop  transfer  func¬ 
tion  data  acquired  at  Test  Point  504  when  Control  Law  N3  was  applied.  The  plots  are  in  polar  form,  allow¬ 
ing  a  simultaneous  comparison  of  magnitude  and  phase  information.  The  numbers  on  the  figure  correspond  to 
frequencies  in  Hertz.  It  is  seen  that  there  is  a  good  agreement  between  the  two  plots  in  both  magnitude 
and  phase.  When  one  considers  that  the  magnitude  and  phase  information  is  influenced  by  mass,  stiffness, 
aerodynamic  and  control  surface  factors,  the  comparison  is  all  the  more  impressive.  The  agreement  appears 
to  deteriorate  as  the  frequency  increases,  but  a  qualitative  resemblance  remains. 

Figures  51,  52  show  a  similar  comparison  for  the  (84  -  83)  response  due  to  the  leading  edge  control 
surface  input.  The  comparison  is  seen  to  be  very  good  for  the  lower  frequency  modes  but  the  higher  modes 
have  a  measured  phase  that  is  90°  ahead  of  the  analytically  predicted  phase.  The  magnitudes  are  seen  to 
compare  quite  well. 

The  preceding  four  figures  dealt  with  responses  to  the  leading  edge  control  surface.  If  responses  to 
the  trailing  edge  control  are  examined,  the  comparison  between  analytical  and  empirical  results  is  not  as 
good.  Figures  53,  54  make  this  comparison  for  (83  -  ai)  under  trailing  edge  excitation  at  Q/Qf  -  1.07. 

For  this  transfer  function,  the  measured  phase  of  the  flutter  mode  is  approximately  45°  ahead  of  the 
analytical  one.  The  higher  modes  differ  considerably  in  magnitude.  The  analytical  mode  at  17  Hz  does  not 
appear  in  the  experimental  result  because  the  higher  modes  have  been  removed  for  the  sake  of  clarity. 

The  comparisons  presented  here  serve  to  explain  ~  number  of  events  that  occurred  during  the  course  of 
the  tunnel  test.  The  most  significant  observation  is  v. at  the  close  agreement  between  the  analytical  and 
the  measured  results  for  the  leading  edge  case  explains  why  the  leading  edge  laws  performed  as  well  as 
they  did.  During  the  test,  the  only  changes  that  were  required  in  the  analytically  derived  control  laws 
were  increases  in  the  overall  gain  levels  by  small  amounts. 

The  empirical  data  also  explains  why  the  N1  and  BAE  control  laws  had  difficulty  with  a  13  Hz  mode. 

The  indication  is  that  the  analytical  phase  information  for  this  mode  was  in  error  to  the  extent  that 
these  control  laws  were  driving  this  mode.  Further  study  is  required  to  determine  why  the  N3  law  did  not 
have  this  difficulty. 

The  primary  conclusion  from  the  trailing  edge  data  is  that  there  are  significant  discrepancies 
between  the  measured  and  analytical  results.  The  least  squares  synthesis  procedure  was  applied  during  the 
test  program  to  modify  the  analytically  derived  N3T  law,  with  the  altered  control  law  showing  some 
improvement. 


Both  the  analytical  and  the  experimental  data  indicate  why  there  has  been  difficulty  in  making  the 
trailing  edge  surface  an  effective  controller.  At  a  dynamic  pressure  which  is  7  percent  above  the  open 
loop  flutter  speed.  Figures  53,  54  show  that  there  is  a  relatively  small  phase  range  in  the  flutter  loop. 
Additional  investigation  shows  this  situation  worsens  as  the  dynamic  pressure  increases,  and  that  there 
may  be  problems  with  higher  modes. 

Detailed  examination  of  the  trailing  edge  data  indicates  that  the  inability  to  accurately  predict 
the  analytical  response,  coupled  with  the  apparently  non-optimal  sensor  selection,  combined  to  make 
trailing  edge  control  difficult.  This  is  unfortunate  since  the  tunnel  tests  indicated  that  this  surface 
is  quite  powerful  in  the  sense  that  it  creates  strong  responses  for  relatively  small  surface  motions. 
However,  based  on  the  experiences  gained  in  the  wind  tunnel  entry,  it  is  now  possible  to  organize  new 
sensor  combinations,  or  to  define  a  new  trailing  edge  surface,  for  the  purpose  of  implementing  trailing 
edge  control  laws  and  achieving  effective  flutter  suppression. 

CONCLUSIONS 

Aircraft  active  wing/store  flutter  suppression  has  drawn  substantial  interest  in  the  technical  com¬ 
munity  dealing  with  aeroelastic  problems.  In  concentrating  on  one  severe  flutter  condition  of  the  fighter 
aircraft  model,  it  was  possible  to  draw  on  the  talents  from  two  continents  to  organize  control  laws  and 
test  them  in  one  facility.  The  differences  in  background  and  flutter  control  philosophy  of  the  contribu¬ 
tors  are  clearly  reflected  in  the  diversified  control  laws  submitted  for  experimentation.  It  is,  there¬ 
fore,  all  the  more  gratifying  that  in  the  final  test,  all  the  control  laws  succeeded  in  serving  their 
intended  functions  as  visualized  by  their  originators. 

The  wind  tunnel  test  demonstrated  a  nunfcer  of  milestone:,  which  are  listed  below: 

-  Applying  a  number  of  control  laws  to  a  severe  flutter  condition  at  M  =  0.80,  the  model  was  tested 
at  Q/Qf  -  1.70,  indicating  a  flutter  speed  improvement  of  approximately  30t.  For  one  control  law, 
a  peak-hold  damping  prediction,  substantiated  by  other  appraisals  of  the  test  data,  indicated  a 
potential  for  the  control  law  to  reach  a  dynamic  pressure  131?,  above  the  open  loop  flutter 
dynamic  pressure. 

-  The  ability  to  switch  from  one  control  law  to  another  at  Q/Qf  =  1.40  was  demonstrated.  This 
feature  is  significant  in  future  adaptive  flutter  control,  where  control  law  changes  at  super¬ 
critical  speeds  are  necessary. 

-  The  ability  to  switch  from  a  leading  edge  control  law  to  a  trailing  edge  control  law,  and  vice 
versa,  was  also  demonstrated. 

-  Test  monitoring  programs  were  developed  to  closely  monitor  the  flutter  control  system  behavior 
in  the  wind  tunnel.  Using  the  programs,  a  control  law  was  monitored  and  improved  by  a  moderate 
amount  at  the  test  site.  Phase  adjustments  were  acconipl  ished  by  changing  the  control  law 
constants,  paving  the  way  for  future  mechanization  of  control  laws  using  digital  computers. 

As  described  in  the  paper,  post  analysis  of  the  test  data  indicated  additional  potential  for  the  flut¬ 
ter  suppression  system,  as  well  as  shortcomings  which  can  be  corrected  in  the  future.  For  instance,  the 
trailing  edge  control  laws  may  be  improved  through  sensor  location  optimization  to  expand  the  phase  mar¬ 
gins.  The  control  law  using  two  surfaces  yielded  experimental  data  which  may  be  applied  to  investigate 
the  unsteady  aerodynamic  intercoupling  effect,  etc.  Additional  processing  of  the  test  data  will  be 
continued  and  reported. 

The  consistent  performance  of  the  active  flutter  suppression  systems  reinforces  our  confidence  in  the 
approach.  Another  technical  step  ahead  involves  the  development  of  adaptive  control  laws  for  flutter  sup¬ 
pression.  The  adaptive  system  is  needed  to  deal  with  the  varying  wing/store  configurations,  or  flight 
conditions,  or  both.  With  the  advancement  in  digital  computers  and  integrated  aircraft  control  system, 
and  through  the  concerted  effort  of  the  technical  community,  it  is  our  vision  that  working  flutter  sup¬ 
pression  systems  will  eventually  be  incorporated  in  modern  combat  aircraft. 
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Figure  1.  YF-17  Wing/Store  Flutter  Suppression  Model,  Configuration  B. 


Figure  2.  Block  Dianram  of  the  Flutter  Detector. 
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Figure  3.  Oscillographs  of  the  YF-17  Model  Transducers  When  a  Moderate  Flutter  Condition  Was 

Encountered.  From  Top  to  Bottom,  the  Six  Data  Traces  Are:  Leading  Edge  Servo  Input, 
Leading  Edge  Surface  Response,  Trailing  Edge  Surface  Response,  Tail  Surface  Root 
Section  8end1ng,  Wing  Root  Section  Torsion,  and  Wing  Root  Section  Bending. 
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Figure  4.  The  Flutter  Stopper  Assembly  in  the  Sparrow  Missile. 
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Figure  5.  Details  of  the  Flutter  Stopper  Assembly. 


Figure  6.  Matched  Velocity  Flutter  Calculations, 

Configuration  B,  M  =  0.8  for  the  Flutter 
Stopper  Inactive  and  Active. 
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Figure  7.  Representative  Block  Diagram 
for  a  Control  Law  Used  During 
the  Tunnel  Test. 


Figure  10.  The  Trailing  Edge  Surface 
Control  Law  (N3T). 
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Figure  8. 
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Figure  9.  The  Three  Transducer  Control  Law  (N3). 
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Leading  Edge  Surface  Control  Law  (Nl). 


F,,F2  represent  the  two  feedback  loops  of  Control 
_  _  Law  (N3)  ,  with  the  K  block  Included. 

Fi,F«,  represent  the  two  feedback  loops  of  Control 
I.aw  (N3T),  with  the  K  block  Included. 

The  Cj  j  elements  are  the  aircraft/servo  trails  ter 
functions,  1  *  transducer  index,  j  ■  servo 
signal  index.  Test  Point  487. 

Figure  11.  The  Two  Surface  Control  Law  (N3P). 
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Figure  12.  Characteristic  Diagram  Results  of 
the  YF-17  Model  With  and  Without 
the  Leading  Edqe  Surface  Control 
Law  Nl,  K  =  0.75. 


Figure  13.  Open  Loop  Nyquist  Plots  of  the  YF- 1 7  Model 
With  Control  Law  (N3),  K  =  0.75. 
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Figure  14.  Open  Loop  Nyquist  Plots  of  the  YF-17  Model 
With  Control  Law  (N3T),  K  =  1.00. 
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Figure  15.  Characteristic  Diagram  Results  of 

the  YF-17  Model  With  the  Two  Surface 
Control  Law  (N3P). 

M  *  0.80,  K.r  ■  0.75,  KTr  -  1.00. 
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Figure  19.  Damping  Trend  of  the  YF-17  Flutter  Suppression  Model  -  Open  Loop 
M  =  0.80. 
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Figure  20.  Damping  Trend  of  the  Northrop  Leading  [dge  Surface  Control 
Law  -  N 1  (TP89-103). 
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Figure  21.  Damping  Trend  of  the  Northrop  Three  Transducer  Control 
Law  -  N3  (TP  117-127). 
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Figure  22.  Damping  Trend  of  the  Modified  Northrop  Trailing  Idge  Surface 
Control  Law  -  N3T  (TP  245-249). 
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Figure  23.  Damping  Trend  of  the  Northrop  Two  Surface  Control  Law  -  N3P  (TP  487-50K). 
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Figure  27.  The  PSD  of  the  Leading  Edge  Surface 
Angular  Displacement  Acquired  at 

Test  Point  104,  Q  =  128  psf. 

Control  Law  Nl. 

Figure  30.  The  PSD  of  the  Wing  Root  Torsion 
Gage  Output  Acquired  at  Test  Point 
64,  Open  Loop,  Q  =  73  psf. 
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Figure  28.  The  PSD  of  the  Leading  Edge  Surface 
Angular  Displacement  Acquired  at 
Test  Point  127,  Q  =  128  psf. 

Control  Law  N3. 
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Figure  31.  The  PSD  of  the  Wing  Root  Torsion 
Gage  Output  Acquired  at  Test  Point 
104,  Q  =  128  psf.  Control  Law  Nl. 


1U"'*  V2 


10  20  30  4  0  SO 

FREQUENCY  (HZ) 

Figure  29.  The  PSD  of  the  Leading  Edge  Surface 
Angular  Displacement  Acquired  at 
Test  Point  211,  Q  =  128  psf. 

Control  Law  BAE,  K  =  1.20. 
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Figure  32.  The  PSD  of  the  Leading  Edge 
Surface  Angular  Displacement 
Acquired  at  Test  Point  610, 

Q  -  128  psf.  Control  Law  N3P. 
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Figure  33.  Transfer  Function  (a3-ai)  Versus 

lSIN,LE-  Test  Point  506,  Q  =  95  psf. 
Control  Law  N3P. 
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Figure  34.  Transfer  Function  (a»-a3)  Versus 
'MN.LE'  Test  Point  506,  0=95  list. 
Control  Law  N3P. 
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Figure  35.  Transfer  Function  ( a 3  -  a t )  Versus 
fIN,TL'  Test  Point  507,  Q  =  95  psf 
Control  Law  N3P. 
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Figure  36.  Transfer  Function  (a* -a,)  Versus 
6 1 n , f p .  Test  Point  507,  Q=95psf 
Control  Law  N3P. 
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Figure  42.  Extracted  Open  Loop  Transfer  Func¬ 
tion  (a3-ai)  Versus  the  Leading 
Edge  Control  Surface  Input  6jn  LE. 
M  =  0.8,  Q  =  95  psf.  ‘ 
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Figure  43.  Extracted  Open  Transfer  Function 
(a^-a,)  Versus  the  Leading  Edge 


Control  Surface  Input  6IN  ,F. 
M  *  0.8,  Q  =  95  psf.  ’ 


Figure  44.  Open  Loop  Nyquist  Plot  for  the 
Example  Presented  in  Figures  39 
Through  43.  The  Test  Condition 
was  M  =  0.8,  Q  =  95  psf,  which 
is  27 %  Above  the  Critical  Open 
Loop  Flutter  Dynamic  Pressure. 


Figure  45.  Open  Loop  Flutter  Results  from  the 
Wind  Tunnel  Test. 
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SUMMARY 

Control  laws  were  calculated,  using  optimal  control  theory,  to  suppress  an  explosive 
wing-store  flutter  case  on  a  YF-17  dynamically  scaled  model.  The  trailing  edge  flap  was 
used  for  flutter  suppression  because  usually  hydraulically  driven  ailerons  are  available 
in  modern  fighters.  The  design  aim  of  1.5  times  the  flutter  dynamic  pressure  could  be 
demonstrated  during  the  wind  tunnel  test.  It  should  be  emphasized  that  no  changes  to  the 
analytically  developped  control  law  were  neccessary  in  the  test  which  proves  that  theory 
has  well  advanced  during  the  last  years. 


INTRODUCTION 

In  addition  to  the  joint  F4  flutter  suppression  program  of  the  U.S.  Air  Force  Flight 
Dynamics  Laboratory  and  the  MBB-Airplane  Division  Control  Laws  were  calculated  by  MBB  for 
active  control  of  an  explosive  wing  store  flutter  case.  These  control  laws  were  imple¬ 
mented  into  a  dynamically  scaled  model  of  the  YF-17  aircraft  by  the  Northrop  Corporation 
Aircraft  Group  and  tested  in  the  16  ft  Freon  Wind  Tunnel  of  the  NASA  Langley  Research 
Center.  Design  aim  was  to  reach  1.5  the  flutter  dynamic  pressure  with  the  control  system 
on.  This  goal  could  be  fulfilled  but  the  necessary  safety  margins  could  not  be  met.  An 
improved  system  was  laid  out  afterwards  analytically  which  would  meet  the  required  mar¬ 
gins.  Since  good  correlation  of  test  results  with  the  first  set  of  control  laws  was  found 
it  is  believed  that  the  improved  system  would  also  behave  as  predicted  in  the  wind  tunnel . 
The  work  performed  by  MBB  during  this  program  was  sponsored  by  the  ZTL-Research  Program 
of  the  German  Ministry  of  Defense.  The  wind  tunnel  program  was  procured  by  the  AFFDL . 

A  complete  description  of  the  whole  program  which  involved  the  cooperation  of  several 
countries  is  given  in  /I/. 


1 .  SYSTEM  DESCRIPTION 

The  flutter  model  used  is  a  30%  scale  half  model  of  the  YF-17  with  an  empty  tip  laun¬ 
cher  rail  and  an  AIM-75  missile  attached  to  a  pylon.  A  picture  of  the  model  is  presented 
in  Fig.  1.  The  model  was  equipped  with  a  hydraulically  driven  leading  edge  and  trailing 
edge  flap  and  three  accelerometers  which  provided  the  feedback  signals  for  our  control 
laws.  The  actuator  transfer  function  is  shown  in  Fig. 2  .  A  ground  resonance  test  was  per¬ 
formed  on  the  model  with  the  tunnel  support  conditions  simulated  121 .  Fig.  3  presents 
the  measured  mode  shapes  of  Mode  I  to  Mode  IV.  The  flutter  calculation  in  Fig.  4  shows 
classical  wing  bending  (Mode  I)  -  wing  torsion  (Mode  II)  flutter.  It  can  also  be  seen 
in  this  figure  that  the  system  reduced  to  a  binary  system  behaves  almost  identical  which 
is  very  useful  for  applying  optimal  control  theory  with  complete  state  vector  feedback. 
All  control  laws  were  laid  out  considering  Mode  I  and  Mode  II  only. 


2 .  ANALYTICAL  APPROACH 

Unsteady  Aerodynamic  Forces 

Unsteady  aerodynamic  forces  were  performed  using  eight  measured  vibration  modes  of  the 
model.  The  rigid  body  modes  were  neglected.  For  the  unsteady  aerodynamic  forces  the  fol¬ 
lowing  simplifications  were  made. 

.  It  was  assumed  that  the  aerodynamic  forces  of  the  tailplane  have  no  influence  on 
the  wing  flutter  mechanism.  Therefore  only  the  wing  aerodynamic  unsteady  forces 
were  introduced  into  the  calculation. 

.  The  lift  of  the  missile  body  and  the  missile-wing  interference  were  neglected. 

The  unsteady  wing  aerodynamic  forces  were  calculated  for  the  eight  measured  ei 
of  the  model  by  the  doublet  lattice  method  / 3 /  and  by  the  lifting  surface  method  <>  ,  . 

The  flutter  calculation  was  performed  at  Ma  =  0.8.  Results  for  the  eight  degrees  <'r 
freedom  system  are  given  in  Fig.  4. 


A 


Because  theory  overpredicts  control  surface  unsteady  airforces  the  trailing  edge 
lift  coefficient  was  multiplied  with  a  factor  of  0.7.  There  are  several  unsteady  pres¬ 
sure  distribution  measurements  in  the  transonic  regime  which  prove  this  correction  /5/, 
/6/. 


Control  Law  Definition 

The  technical  approach  to  find  the  control  law  is  described  in  /7/. 

To  get  the  optimal  control  law  k0pt  the  quadratic  performance  criterion  is  minimized 

I  [(|.|'  [O]  |.|  •  ..K  .  )  a. 

0 

Q  is  a  weighting  matrix  found  by  trial  and  error  by  using  a  screen  together  with  the 
computer.  R  is  a  scalar,  to  be  selected,  because  there  is  only  one  control  surface.  X  is 
the  state  vector  and  Xi  is  the  actuator  input  signal. 

Minimizing  (J)  leads  to  the  optimal  control  laws 


with 

where  P  is  the  steady  state  solution  of  the  Matrix  Riccati  equation 

H-HM-M'H 

Because  there  is  no  direct  access  to  the  state  vector  x,  linear  combinations  of  the 
sensor  outputs  can  be  used  to  get  the  state  terms 

fyi  -  rc](x| 

where  y  is  the  sensor  output  and  C  a  transformation  matrix. 

The  following  sketch  shows  the  block  diagram  for  this  model. 


(x.! 

(*l 

1  u 

B 

r 

- V 

J 

1 

r 

| 

A 

I 

I  AIRPIANF.  DYNAMICS 

( _ 

K 

■;  AM  WITH 

!*! 


KM 


'ST 


lYl 


The  control  vector  Km  related  to  the  measured  vector  y  can  be  defined  by  the  inverse 
transformation  matrix  C-1 


K,r=  s  % y  c  c-'j 


No  problem  arises  to  perform  equation  aix>ve  if  the  complete  state  vector  is  fed  back  thus 
leading  to  a  quadratic  form  of  tne  transformation  matrix  C. 

The  flutter  case  can  be  described  by  two  modes 


.  first  wing  bendlnq  g( 

.  first  wing  torsion  +  stoic  pitch  q. 

These  modes  can  be  measured  by  three  accelerometers  on  the  wing.  For  the  realization 
of  the  control  law  it  was  found  that  contribution  of  the  states  related  to  the  aileron  is 
small  and  therefore  no  feed  back  of  those  states  is  necessary. 

The  block  diagram  of  the  control  law  realized  for  tunnel  testing  is  shown  in  Fig.  r> . 
The  factors  ka,  Kft  are  related  to  t  lie  wing  torsion  and  K<t‘  ,  K  ^  arc  related  to  the  wing 
bendinq.  It  was  also  possible  to  change  the  total  gain  and  the  total  phase  of  the  actu¬ 
ator  input  signal  Xj. 


The  control  law  was  found  by  using  only  two  generalized  coordinates  (wing  bending 
wing  torsion)  as  feedback.  When  the  total  vibration  system  is  used  then  there  is  always 
the  possibility  that  adverse  signals  are  picked  up  because  the  measured  vector  contains 
contributions  from  all  modes.  A  flutter  calculation  plot  is  presented  in  Fig.  6  which 
shows  that  Mode  IV  (fuselage  first  bending.  Fig.  3)  becomes  unstable.  There  is  consider¬ 
able  motion  of  the  outer  wing  on  this  mode.  The  Nyquist  plot  of  Fig.  7  shows  a  perfect 
circle  for  the  torsion  mode  with  sufficient  gain  and  phase  margin  (6  db  *_  60°)  .  It  also 
shows  a  mode  at  about  12.8tlz(Mode  IV)  to  be  unstable  because  the  point  -1  is  encircled 
from  the  left  side.  The  bandpass  filter  of  Fig.  8  was  incorporated.  The  effect  of  this 
filter  was  compensated  only  at  6  Hz  by  shifting  the  control  law  forward  about  80°  (which 
means  just  other  factors  for  Kij>,  Ki}>,  Ke,  Ke  of  Fig.  5)  . 

A  Nyquist  plot  for  this  new  control  law  is  presented  in  Fig.  9  which  shows  Mode  IV 
to  be  stable  but  the  circle  at  6  Hz  is  deteriorated  (there  is  less  phase  margin  now!). 

An  adverse  behaviour  of  this  control  law  is  also  shown  in  the  v-g  plot  of  Fig.  10.  The 
flutter  speed  is  reduced  about  10  kts  (compared  to  Fig.  6)  and  the  damping  behaviour  is 
much  more  explosive  because  now  the.  bending  mode  is  fluttering. 


3.  WIND  TUNNEL  TEST  RESULTS 

The  peak  hold  spectrum  method  was  extensively  used  to  find  the  model  flutter  speeds 
with  and  without  flutter  control.  The  method  provides  power  spectrum  density  values  for 
different  signals.  The  input  signal  is  wind  tunnel  noise. 

Typical  examples  of  this  method  are  plotted  in  Fig.  11.  The  inverse  of  the  amplitude 
for  the  two  peaks  of  Fig.  11  is  plotted  in  Fig.  12  versus  speed.  The  model  flutters  at 
the  speed  where  this  value  is  zero.  Fig.  11  indicates  the  same  behaviour  as  shown  in 
the  v-g  plot  of  Fig.  10.  Whereas  the  torsional  mode  is  suppressed  with  increasing  speed 
the  bending  mode  becomes  unstable.  Fig.  13  shows  the  model  flutter  speed  versus  total 
gain  of  the  flutter  suppression  system.  From  this  figure  it  can  be  deduced  that  an  in¬ 
crease  in  gain  above  0.7  increases  the  flutter  speed  very  little  and  would  only  make  the 
systems  very  sensitive  to  disturbances.  The  phase  margins  of  +  60°  could  not  be  fulfilled 
with  this  control  law  whereas  the  6  db  gain  margin  can  almost  be  met. 

Fig.  14  shows  good  correlation  of  test  results  with  analytical  predictions. 


4.  IMPROVED  FLUTTER  SUPPRESSION  SYSTEM 

After  consideration  of  wind  tunnel  test  results  and  a  review  of  control  law  finding 
procedures,  two  major  deficiencies  were  found 

.  the  band  pass  filter  used  had  large  phase  shifts  in  the  range  from  4  Hz  to  8  Hz 
(about  90°  phase  lag) 

.  it  was  not  included  in  the  optimization  procedure  but  just  compensated  for  one 
distinct  frequency. 

A  band  pass  filter  was  designed  which  had  almost  the  same  attenuation  behaviour  as 
the  one  used  in  the  test,  but  had  much  less  phase  shift  in  the  range  4  Hz  to  8  Hz  (about 
20°  phase  lag).  The  filter  behaviour  is  shown  in  Fig.  19. 

The  filter  transfer  function  was  multiplied  with  the  actuator  transfer  function  and 
approximated  at  the  design  frequency.  So  the  filter  function  was  already  included  when 
optimal  control  theory  in  the  time  domain  was  applied  and  new  control  laws  were  found. 


Analytical  Results  for  the  improved  System 

Having  found  the  new  control  vector,  a  Nyquist  diagram  was  made  at  the  design  speed 
of  280  kts.  This  diagram  (Fig.  16)  shows  now  a  perfect  circle.  The  v-g  plot  of  Fig.  17 
also  shows  a  much  less  critical  damping  behaviour  (compared  to  Fig.  10).  It  is  always 
the  torsion  mode  which  becomes  unstable. 


Figures  18,  19,  20  show  the  results  of  sensitivity  studies  were  the  system  total 
gain  and  the  phase  shift  of  the  actuator  input  signal  were  investigated.  Again  it  is 
proven  that  the  best  control  law  is  the  one  produced  by  applying  optimal  control  theory 
(K  =  1.0,  cp  «  0)  . 

Fig.  18  shows  a  stable  system  with  a  flutter  speed  of  390  kts  for  K  =  1.9.  with  a 
6  db  qain  margin  300  kts  flutter  speed  could  be  achieved  which  is  above  our  design  goal. 
The  phase  margins  are  still  below  the  desired  _*  60°  phase  shifts  as  is  presented  in 
Fig.  19  and  Fig.  20.  For  an  operational  system  one  would  have  to  control  the  phase  of 
the  total  system  very  closely  -  ixissibly  bymeans  of  self  measuring  and  self  adjusting. 
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1G.  1  YF-17  MODEL  AND  MODEL  SUSPENSION 
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FIG.  2  TRAILING  EDGE  FLAP  ACTUATOR 
TRANSFER  FUNCTION 


,i  ;■  :l:  nrt 

UtWT 

|t 

- 

i ;  V 

M 

'si  TORSION 

;  \r  • 

^ - — - 

BINDING  W  ' 

'f 

fM. 

\  | 

I 

>M 

\  ' 

IN 

\  " 

\ 

—  -<MIC  1  \ 

i  -l 


FIG.  3  MEASURED  MODE  SHAPES  OF  THE  MODEL 
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FIG.  A  FLUTTER  CALCULATION  WITH  EXTENDED  AND  REDUCED  SYSTEM 
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G.  5  BLOCK  DIAGRAM  OF  THE  CONTROL  LAW 
TESTED  IN  THE  WIND  TUNNEL 
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FIG.  6  V-G  PLOT  OT  CONTROL  LAW  WITHOUT  BANDPASS 
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FIG.  13  FLUTTER  SPEED  VS  GAIN 


FIG.  1A  COMPARISON  OF  TEST  RESULTS 
AND  PREDICTIONS 
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RESUME 

Une  optimisation  d'une  loi  de  controle  de  flottement  a  ete  recherche*-  pour  1*  »  an  d*  . «  ma-p.itt* 
de  1 'avion  YF17  £quipe  d'une  charge  en  bout  d'alle  (configuration  H). 

La  loi  de  controle  appliquee  est  basee  sur  It-  princlpe  de  la  reinjection  d'un*  tlgt  it*  s>.t 
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qul  pouvalent  exister  entre  les  modes  et  les  frequences  utilises  dans  le  >  alnil  ave*  1  »•*.  moil*-*,  it 

frequences  reels  lors  du  passage  en  soufflerie. 
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SUMMARY 

A  f  lutt  er  cent  rol  law  was  calm  b*  ed  f  or  a  YF  1  7  model  hjh  tped  with  *  -t  t  «  at  i  ■  «  •  i  • 

The  control  law  used  stiffness  Injection  on  the  wing  by  the  helj  !  the  unst  a.  r 

Induced  by  a  control  surface. 

In  the  present  case  Just  the  leading  edge  control  stir  fat*  and  only  ru  a 

the  wing  close  to  the  nodal  line  of  the  {'Itching  mode  of  the  external  stot*  ,  writ  ,s.,‘ 

During  the  wind  tunnel  tests  the  nominal  control  law  was  used  without  -vi:,  ,  •  ,  .  *  ; 

differences  between  eigen  modes  and  frequencies  Introduced  in  the  tabulations  uii«;  >  :*•;  •«; 

c  les  existing  on  the  model  itself  when  mount  «d  In  the  wind  t  untie  1. 

The  control  at  constant  Mach  number  M  ■  ;  *  i  m  1 1  t  <*!  t  <  in-  r  *  a  si  s  »•  «  ’  ar  •  w»  .  .  ■  .  .  •  .»r> 

cal  pressure  and  als<<  to  reduce  considerably  the  exploslvitv  of  the  I  hittrt  p».»  •.  -rr-.-i. 


Comparison  between  theory  and  tests  is  correct  qualitatively  an!  quant  Hat  1  -  • 


INTRODUCTION 

Dans  le  cadre  d'un  travail  de  cooperation  entre  U.S.  Air  Force  et  3  pays  europeens  (Grande- 
Bretagne,  Allemagne  Federale  et  France)  il  a  ete  decide  que  chaque  pays  experlmenterait  une  loi  de  controle 
de  flottement  calculee  par  lui-meme  et  appliquee  a  une  maquette  dynamiquement  semblable  de  1 'avionyF17 .  Une 
charge  en  bout  d'aile  etait  lestee  de  faijon  fictive  et  permettait  d 'obtenir  un  flottement  tres  violent  a 
Mach  0,80  dans  la  soufflerie  16  pieds  au  freon  de  la  NASA  a  Langley  Field. 

La  maquette  et  les  essais  etaient  sous  la  responsabilite  de  la  firme  Northrop  (Fig.  1). 

Du  cote  fran^ais  l'etude  s'est  deroulee  en  trois  phases  : 

a)  Montrer  que  l*on  etait  capable  de  prevoir  correctement  le  comportement  aerodynaraique  de  la 
maquette  jusqu'au  flottement  sans  controle. 

b)  Determiner  une  loi  de  controle  jugee  valable  pour  ce  type  de  flottement  et  la  comparer  avec 
les  resuitats  experimentaux. 

c)  Calculer  une  loi  optimisee  a  partir  des  donnees  aerodynamiques  determinees  en  a)  et  b) . 
a)  Ca leu 1  de  flottement 

Les  calculs  de  flottement  ont  ete  effectues  a  Mach  0,80,  qui  etait  le  Mach  choisi  pour  les  essais. 

Les  hypotheses  suivantes  ont  ete  introduites  : 

.  Les  six  premiers  modes  de  la  maquette  ont  ete  consideres  (2  modes  "rigides"  de  suspension, 
translation  et  tangage  et  4  modes  de  deformation).  La  gamme  de  frequence  etudiee  s'etendait  de  0  a  15  Hz. 

.  Les  t <*r.  tb  aerodynamiques  instationnaires  ont  ete  calculees  par  une  methode  de  doublets  subso- 
niques,  les  surfaces  considerees  etant  l'aile  et  l'empennage  horizontal.  Aucune  force  n'a  ete  introduite 
sur  le  pylone  ou  l'engin  AlM75,de  la  meme  fa^on  les  terraes  d *  interaction  volume/empennage  et  engin- 
mat/voilurc  ont  ete  negliges. 

.  Le  couplage  est  pr Inc ipalement  du  aux  modes  de  flexion  voilure  et  tangage  engin  qui  induit  une 
torsion  sur  l'aile. 

.  la  figure  2  donne  devolution  des  parametres  de  frequence  et  d 'amor tissement  des  differents 
modes  en  function  de  la  densite  du  fluide  dans  la  soufflerie  au  freon  a  Mach  0,80.  Il  faut  noter  la  tres 
forte  explosivite  du  phenoroene,  e'est-a-dire  la  pente  du  mode  critique  (ici  le  mode  de  tangage). 

dp 

.  Les  essais  comparatifs  qui  sont  portes  sur  la  figure  3  montrent  que,  malgrt*  quelques  differences 
dans  les  donnees  entre  calculs  et  essais,  1 'allure  du  phenomene  et  les  ordres  de  grandeur  sont  bien  prevus. 

.  Les  amor t Issements  et  frequences  ont  ete  determines  exper imentalement  a  partir  de  1 ' enregistre- 
ment  des  couches de  densite  spectrale  moyennees  durant  30  secondes  en  travail lant  avec  des  spectres  redon- 
dants  des  reponses  a  la  turbulence  naturelle  de  la  soufflerie. 

.  Differents  capteurs  privilegiant  soit  le  mode  de  flexion,  soit  le  mode  de  tangage  de  l'engin, 
ont  ete  utilises  pour  suivre  1 'evolution  de  ces  parametres. 

Cette  premiere  partie  montrait  qu’en  ce  qui  concerne  les  forces  aerodynamiques  instationnaires 
sur  l'aile,  les  hypotheses  faites  etaient  valables,  tout  au  moins  au  nombre  de  Mach  M  =  0,80. 

h.M  Lalcul  de  la_Joi_de_  controle 

Les  equations  du  flottement  s'ecrivent  en  notation  matricielle  : 

{Z(l<)]^fM  =  "Cc ( U )]  ® 

od  Z(k,  est  la  mat  rice  d' impedance  du  system*  qui  s'eorit  : 

[*<„,]  =  {-[>]-’  +  W+  % 

dans  une  base  de  modes  propres,  |X  et  ^  ( respec t i vement  matrice  des  masses  et  des  rigidites  general isees) 
sont  d  iagonales  .  )  f  A(  ^  )  et  B(kJ  sont  les  parties  reelles  et  imaginalres  des  forces  aerodynamiques 

Instationnaires  dues  a  un  mouvement  de  I  *a i  1  eron  (  C(  j<  ))  on  a  un  mouvement  de  l'aile^Af^  )  et  B  d< )  ] 

1  *n  t  i r «■  :  q  =  -[Z‘'HC]  ® 


,  t  [w][q]=  -K][Z]  lC  ]  6 


^  Af  ♦■■it  ia  mat  rift-  des  modes  propres  sans  vent  au  1/4  uvant  !*•  la  tranche  d'aile-  <n  In*  d* 


En  posant  Z 


:  i  =_|\Vl'IZ(k,f|C(k)UTts) 


C$L)  est  la  colonne  des  forces  induites  sur  1 'aile  par  la  surface  de  controle,  la  partie  aerodynamique 
de  ces  forces  s’ecrit  : 

S,  X,  |  o  o  o 

0  0  md 

0  3  I  1  _  L  0  0  nd J 


s. 

y. 

_w_ 

0 

0  ii 

C, 

C„ 


La  partie  relative  aux  forces  d'inertie  de  la  gouverne  a  etc  negligee  apres  un  calcul  prelimi- 
naire  comparatif. 

l-a  fonction  ^(5)  a  ete  calculee  pour  200  valeurs  de  frequence  dans  la  gamine  0,20  Hz.  Les  forces 
aerodynamiques  des  fonctions  A,  B,  C  ont  ete  interpolees  pour  chaque  frequence  reduite  a  partir  du  calcul 
exact  effectue  pour  7  valeurs  discretes  de  k  (L* interpolation  est  du  type  parabolique)  : 


0  =  K»  [  S'* . ] 


«1  sont  les  amplitudes  des  deplacements  du  point  de  la  tranche  d  'aile  comportant  la  gouverne  de 
controle  et  situe  au  quart  avant  de  eelle-ci.  Les  deplacements  etant  consideres  dans  les  six  modes  intro- 
dults  dans  le  calcul  (sans  vent). 

^tisont  les  deplacements  du  point  de  controle  dans  le  vent  (point  86,  Fig.  1). 


Fn  boucle  ouverte  on  peut  ecrire  pour  le  systeme  avec  controle  : 

"*O.L  ~  “  Coa  Ts  Tc 

est  la  fonction  de  transfert  de  la  loi  de  controle  etT^^  — 


_<^sTc 


A  —  cov  ‘s  TC 


en  boucle  fermee. 


l.a  fonction  T q  s'ecrlt  : 

Tc  ±  G  T,(g.  +  Gj  T4')T(,  radian/m/s2 
N  *  f litre  de  bands  pour  limiter  le  nombre  de  modes  a  controler  : 

T‘  1  1 


0.07958S 


1  +  0, 8 105S  +  0,0001  75<1S**  I  +  0,0245S  +  0,0001759S^  1  +  0,07  9  58S 

est  un  pseudo- integrateur  qui  coupe  les  frequences  basses  : 

0,282S 

•  - - 

1  ♦  0,282S  +  i',212S‘ 

T  5  est  la  fonction  de  transfert  du  verin  et  do  la  scrvi-valvi  f-'urr.is  par  I’SAF 

G j  =  0,1  C.2  =  Gj  *  1 


Le  schema  bloc  de  la  figure  A  montre  le  prineipc  de  la  loi  do  controle.  Far  ajustement  des  gains 
Gj  et  IF,  la  phase  globale  de  la  fonction  de  transfert  T<~  est  eg  a  l  e  a  180°  a  la  1  -equence  critique  de  *'  H? 
du  mode  a  controler.  La  reinjection  des  forces  sur  l'aile  est  alors  une  rigidity  pure  it  negative,  compte 
tenu  des  hypotheses  sulvantes  introduites  dans  le  calcul  : 

-  les  forces  aerodynamiques  induites  par  la  gouverne  de  controle  n'agissent  de  fat, on  senslMi  qtn-  dans  la 
tranche  d'aile  ou  est  situee  cette  gouverne  ; 

-  les  parties  replies  des  forces  aerodynamiques  creees  par  cette  gouverne  dependent  peu  de  la  frequence 

reduite  k  s.  dans  un  intervalle  de  quelques  Hertz  (de  5  a  10  Hertz  par  cwniplv) 

-  les  parties  Imaginaires  des  forces  aerodynamiques  sont  falbles  devant  les  parties  reel les  en  sorte  qtie  le 
dephasage  qu'elles  introdulsent  peut  etre  neglige. 

i:,'B  ‘1  Iff,  rentes  hypotheses  ont  ete  verifiers  par  des  essals  en  soufflerte  dans  lesquels  on  me  sura  it 
les  Champs  de  press  Ion  instat ionnaire  crees  par  one  nouvernc  exciter  harmoniquement  a  d  1  f  f  .-rentes  frequences. 

t  (I. a  forme  en  plan  des  maquettes  qui  ont  servi  a  ces  mesures  elait  voisine  de  relle  de  l'aile  de 

1  av ion  yF  1 7 ) . 

11  re  so  1 1  e  de  ces  hypotheses  quo  les  forces  Induites  par  la  gouverne  sont  prat iqucm.  nt  en  phase 
aver  le  mouvement  de  rello-ci  et  done  (su leant  le  slum-  cholsl)  en  phase  o„  en  opposition  aver  le  mauve 
moot  du  point  do  control**. 


n’oi'j  I’etfet  do  rigid  l  tv  nvg.it  ivo  qui  a  t*u>  utilise 
tli-xlon  tn  rvduisant  l' importance  des  termes  de  couplagc  par 
d-mc  .1  reduire  l  vs  rlsquvs  dv  flottement. 


<t  qui  t  vmi  .1  abaissvr  la  frequence  du  modi  d« 
separation  des  frequences  critiques  et  tend 


A 


b.2)  Comparalson  des  resultats  d'essais  avec  la  theorie  -  Cas  avec  con t role 


Les  valeurs  exper imentales  des  amortissements  et  des  frequences  des  modes  intervenant  dans  le 
phenomene  de  flottement  ont  ete  obtenues  experimentalement  en  analysant  les  courbes  de  densite  spectrale 
pour  differents  capteurs. 

Une  tentative  d ' interpretation  faite  avec  la  "peak  Hold  Method"  (enveloppe  des  valeurs  maximales 
pour  chaque  frequence  des  courbes  de  densite  spectrale)  a  fait  apparaltre  une  large  dispersion  des  araortis- 
sements.  Ceci  peut  etre  du  aux  raisons  suivantes  : 

-  les  valeurs  maximales  au  cours  du  temps  de  la  courbe  de  densite  spectrale  dependent  de  deux  paramdtres 
d'entree  : 

-  du  niveau  de  turbulence  de  la  soufflerie,  qui  est  fonction  de  la  pression  generatrice  et  du 
n ombre  de  Mach  ; 

-  de  1 ' amor tisseraent  global  du  systeme  (aerodynamique  +  structural)  qui  est  le  seul  aui  doit 
etre  considers  ; 

Le  premier  parametre  etant  ni  constant,  ni  connu,  fausse  la  mesure,  d'autant  plus  que  l'on  veut 
proceder  a  une  extrapolation  sur  de  faibles  valeurs  d 'amort issemen t  pour  obtenir  I#1  point  dp 
f lot  tement . 

-  entin,  la  tonction  inverse  du  maximum  de  la  courbe  de  densite  spectrale  est  une  fonction  inconnue,  non 
lineaire  en  general,  qu'il  est  difficile  d ' extrapol er . 

La  figure  5  presente  le  calcul  avec  controle,  elle  doit  etre  comparee  a  la  figure  2  (calcul 
sans  controle). 

La  comparaison  des  resultats  d'essais  et  de  calcul  avec  controle  est  portee  sur  la  figure  6. 

On  peut  noter  1 'accord  correct  qui  e  xiste  quant  a  l'allure  des  courbes  theoriques  et  experiment 

tales . 

L'explositive  du  phenomene  est  la  meme  dans  les  deux  cas  et  a  ete  cons iderablement  reduite  par 
rapport  au  cas  nominal  sans  controle.  La  limite  de  la  pression  dynamique  que  1  'on  pourrait  atteindre  en 
augraentant  encore  la  pression  dynamique  ne  presente  pas  un  grand  interet  si  on  considere  que  la  courbe 
d 'amortissement  du  mode  critique  est  maintenant  tres  plate.  On  rencontrerait  une  entree  en  flottement  tres 
progressive  et  done  mal  definie. 

II  faut  remarquer  sur  la  courbe  des  frequences  que  l'effet  recherche  d '  introduct  ion  d 'une  rigid  itt- 
negative  sur  le  mode  de  flexion  est  similaire  en  calcul  et  en  essais. 

L'amortissement  de  structure  a  ete  soustratt  des  courbes  d 'amort issemen ts  presentees,  v  les  ne 
comportent  done  que  la  partie  d ' aerodynamique  instat ionnaire. 

c )  Amelioration  de  la  lot  de  controle 


A  partir  de  ces  resultats  on  peut  avoir  une  confiance  raisonnable  dans  le  calcul  des  forces  aerody- 
namiques  agissant  sur  l'aile  et  la  gouverne  de  controle.  Une  optimisation  de  la  loi  de  controle  peut  alors 
etre  appliquee. 

L 'amelioration  a  ete  obtenue  par  un  calcul  a  posteriori  en  reduisant  la  variation  de  phase  en 
fonction  de  la  frequence  dans  la  fonction  relative  a  la  loi  de  controle. 

Le  calcul  a  ete  effectue  a  partie  d'une  loi  ideale  dont  le  module  serait  constant  et  egal  a  un 
radian  de  rotation  de  la  gouverne  pour  1  metre  de  deplacement  du  point  de  controle  et  en  prenant  d 'autre 
part  des  pentes  de  variations  de  phase  en  fonction  de  la  frequence  egales  a  : 


1) 

bf 

>  0°/Hz 

2) 

h'? 

-  -7,5” 

3) 

-  - 1 5°/H 

Dans 

tous  les 

cas  1)  2 

phase  constamment  nulle  ou  egale  a  180° 


4) 


at 

hr 


=  -20°/Hz 


Dans  tous  les  cas  1)  2)  3)  et  4)  la  phase  est  calee  egale  a  180®  pour  la  frequence  critique  de 

b  Hert2. 

Les  courbes  devolution  de  frequences  et  d  'amor  t  issemen  t  pour  les  differents  cas  de  controle  sont 
donnees  pour  les  modes  de  flexion  (mode  n°  2)  et  de  tangage  engin  (mode  n°  3).  Bien  entendu  le  calcul  complet 
A  six  modes  a  et£  effectue. 

_n  . 

On  volt  que  pour  w  — w  (phase  nulle  dans  toute  la  bande  de  frequence)  il  n'existe  plus  de 
possibilite  de  flottement,  les  modes  de  structures  ne  sont  jamais  couples,  les  resultats  de  calcul  sont  iden- 
tiques  a  ceux  que  l'on  pourrait  obtenir  en  faisant  le  calcul  mode  par  mode  (calcul  a  1  degre  de  1 iberte) . 

Lorsque  la  pente  de  phase  atteint  -20®/Hz  on  se  retrouve  pratiquement  dans  les  conditions  dc  l'essal 
et  le  flottement  reapparalt  avec  un  gain  sur  la  pression  dynamique  dp  602  du  au  controle. 
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CONCLUSION 


On  a  pu,  dans  cette  etude,  valider  dans  une  premiere  phase  les  hypotheses  restrictives  dans  le 
calcul  des  forces  aerodynamiques  instationnaires  aussi  bien  en  ce  qui  concerne  les  termes  relatifs  a  l'alle 
que  ceux  afferents  a  la  surface  de  controle. 

La  loi  utilisee  pour  reduire  le  flottement,  consistant  a  introduire  un  terme  de  rigidite  sur  1 'un 
des  modes  de  la  structure,  s'est  averee  correcte.  La  comparaison  theorie-experience  est  en  bonne  concor¬ 
dance. 

II  est  ensuite  facile,  partant  d  'une  fonction  de  transfert  connue  pour  la  loi  de  controle,  de 
l'ameliorer  par  le  calcul  en  reduisant  la  pente  de  la  courbe  de  phase  fonction  de  la  frequence. 

Toutefois,  cette  methode  ne  s 'attache  pas  a  obtenir  une  optimisation  sur  les  marges  de  gain  ou  de 
phase  du  controle  en  fonction  de  parametres  comme  la  pression  dynamique  ou  le  nombre  de  Mach. 

Une  loi  de  controle  numerique  pilotee  par  ces  memes  parametres  s’averera  necessaire  dans  le  futur 
pour  des  applications  portant  sur  des  avions. 
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Fig.  /  -  Ft  7  maquette.  Configuration  B. 


Fig.  2  FJ  7  maquette  de  flottement.  Thtorie  sans  contrfile. 
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SURFACE  OE 
CONTR6LE  DU 
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O  (  R5  F) 
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Emois  av*c  control* 


M7 

-a»0  V. 123mA 


O 
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Th«pri#  ovrc  control 


Fig.  6  Comparator!  essgis-iohuh  urn  hi  de  <  ontrdie 
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WIND  TUNNEL  TESTS  ON  A  FIGHTER  AIRCRAFT 
WING/STORE  FLUTTER  SUPPRESSION  SYSTEM 
-  THE  B.Ae  CONTROL  LAW 
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SUMMARY 


A  B.Ae  control  law  designed  using  analytical  data  was  tested  on  a  model  of 
the  YF-17  in  the  NASA  Langley  16'  wind  tunnel  and  succeeded  in  meeting  the  require¬ 
ment  to  increase  the  flutter  dynamic  pressure  by  70%  at  M  =  0.8. 

Because  the  store  configuration  chosen  for  these  tests  produced  a  flutter 
mechanism  which  was  both  hard  and  near  to  classical  frequency  coalescence,  it 
was  necessary  to  use  a  multiple  output  control  law  to  achieve  the  stability 
margin  objectives  of  T60°  and  +6dB. 

The  B.Ae  control  law  was  designed  using  a  novel  procedure  which  provides 
these  stability  margins,  uses  minimum  control  surface  movement  in  turbulence  and 
can  be  designed  using  either  analytical  or  empirical  data.  Two  wing  tip  acceler¬ 
ometers  and  a  leading  edge  control  surface  were  used. 

Empirical  open  loop  transfer  functions  obtained  during  the  test  showed 
that  the  analytical  data  overestimated  the  response  of  the  flutter  mode  to  leading 
edge  control  surface  excitation. 

INTRODUCTION 


The  United  States  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  sponsored 
the  design,  construction  and  wind  tunnel  testing  of  a  YF-17  model  to  demonstrate 
active  control  of  store  flutter. 

Northrop  Corporation  were  contracted  to  design  the  model  and  the  active 
flutter  control  laws,  and  to  test  these  in  the  16'  Transonic  Dynamics  Tunnel  of 
NASA  Langley  Research  Centre. 


After  an  initial  phase  of  wind  tunnel  testing  (ref.  1),  information  exchange 
agreements  were  drawn  up  between  the  AFFDL  and  each  of  three  European  countries 
(U.K.,  Germany  and  France)  whereby  active  flutter  control  laws  would  be  designed 
in  each  of  these  three  countries  for  the  YF-17  model  to  be  tested  along  with 
control  laws  designed  by  Northrop. 

These  control  laws  were  tested  in  October  1979  (ref.  2)  and  this  report 
briefly  describes  the  work  done  by  British  Aerospace.  It  is  in  three  parts: 
the  basic  properties  of  the  model  as  shown  by  analytical  studies,  control  law 
design  and  the  results  of  the  wind  tunnel  testing. 

WIND  TUNNEL  MODEL 

Fig.  1  shows  the  30%  scale  wind  tunnel  half  model.  It  is  supported  vert¬ 
ically  by  soft  springs  and  attached  to  the  tunnel  wall  by  roll  bars  so  that  only 
symmetric  motion  of  the  aircraft  is  represented.  The  wing  and  fuselage  structuie 
is  dynamically  representative  of  the  full  scale  aircraft. 
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Two  control  surfaces  (one  leading  edge  and  one  trailing  edge)  aie  available 
for  active  flutter  control.  Each  is  moved  independently  by  its  own  hydraulic 
actuator. 

For  the  October  1979  tests,  there  were  four  accelerometers  on  the  outer 
wing  available  for  active  flutter  control.  Also  the  model  had  a  pylon  mounted 
AIM-7S  missile  and  an  empty  wing  tip  launcher  rail  (i.e.  store  configuration  B). 

3.  DESIGN  OBJECTIVE 

The  design  objective  laid  down  by  the  AFFDL  was  to  increase  the  flutter 

dynamic  pressure  by  70%  with  stability  margins  oi  i no  and  TodR  using  the  existing 

control  surfaces  and  accelerometers. 

4.  BASIC  PROPERTIES  OF  THE  MODEL 

The  analytical  studies  on  the  wind  tunnel  model  were  made  using  two  sets 
of  aeroelastic  data  provided  by  the  AFFDL  and  Northrop  Corporation.  The  first 
(a  10  degrees  of  freedom  system)  was  based  upon  seven  calculated  normal  modes 
of  the  model.  The  second  (a  13  degrees  of  freedom  system)  was  produced  later 
and  was  based  upon  ten  GVT  modes  modified  to  give  agreement  with  the  critical 
flutter  dynamic  pressure  (QF)  measured  in  earlier  wind  tunnel  tests  (ref.  1). 

4. 1  Flutter  Curves 

Fig.  2  shows  the  basic  flutter  curves  for  the  10  d,  of  f.  system.  The 
flutter  root  is  that  associated  with  the  air  off  fundamental  wing  torsion  mode. 
Below  the  flutter  dynamic  pressure  (QF)  the  flutter  root  damping  does  not  exceed 
1%  critical.  However,  above  QF  the  damping  falls  rapidly  to  about  -10%  critical 
at  1.7QF  i.e.  the  flutter  is  hard. 


Fig.  3  shows  the  flutter  curves  for  the  binary  system  consisting  of  the 
fundamental  bending  and  torsion  modes.  The  frequency  and  damping  curves  are 
identical  to  those  of  the  corresponding  roots  in  fig.  2  with  an  almost  identical 
QF.  This  shows  that  the  flutter  mechanism  is  pure  binary. 

Also  shown  in  fig.  3  is  the  effect  of  ignoring  the  out  of  phase  aerodynamic 
generalised  forces  (  [&]  »  [o]  ),  This  shows  the  classical  frequency  coalescence 

flutter  (ref.  3  and  4)  with  the  flutter  speed  not  very  much  changed. 

Therefore  fig.  2  and  3  show  that  we  are  attempting  to  perform  active  flutter 
control  of  hard,  near  frequency  coalescence,  binary  flutter. 

Open  Loop  Transfer  Functions 


4.2.1  Input  to  the  Leading  Edge  Control  Surface  Actuator 

Fig.  4  shows  the  open  loop  transfer  function  (o.l.t.f.)  at  1 . 7QF  for  input 
to  the  leading  edge  control  surface  actuator  and  velocity  output  from  the  outboard 
aft  (0/A)  accelerometer  position  computed  using  the  10  d.  of  f.  system  data. 

There  are  two  significant  features  which  were  seen  from  the  o.l.t.f.  for 
outputs  at  lOO  positions  over  the  wing  surface.  These  werei  all  anticlockwise 
lobes  are  narrow  and  have  the  same  mean  phase  angle,  and  the  wing  tip  transducers 
sense  the  least  response  in  the  modes  not  involved  in  the  binary  flutter. 

This  means  that  a  simple  control  law  using  a  single  transducer  will  have 
inadequate  phase  margins,  and  a  control  law  which  uses  the  two  outboard  transducers 
will  have  least  problems  avoiding  closed  loop  instability. 

4.2.2  Explanation  of  Open  Loop  Transfer  Function  Properties 

Fig.  5  and  ref.  5  show  that  if  the  dynamic  system  is  simplified  toi 

a  binary  flutter  system  without  out  of  phase  aerodynamics 

excited  by  only  the  in  phase  aerodynamics  of  a  control  surface, 

the  o.l.t.f.  between  the  control  surface  angle  (3  and  the  displacement  anywhere 
on  the  structure  W  is  the  product  of  two  transfer  functions  each  of  which  has 
real  values  for  all  frequencies. 

As  a  result,  the  velocity  o.l.t.f.  w/(S  will  have  values  only  on  the 
imaginary  axis,  starting  and  finishing  at  the  origin.  This  means  that  a  single 
input/single  output  control  law  for  such  a  system  will  have  zero  phase  margins. 
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Fig.  6  shows  what  happens  when  the  YK-17  model  is  represented  in  this 
simplified  manner.  The  velocity  output  at  O/A  due  to  the  leading  edge  control 
surface  motion  gives  an  o.l.t.f.  which  travels  up  then  down  the  positive  imaginary 
axis  as  frequency  increases.  However,  when  the  trailing  edge  control  is  used,  the 
o.l.t.f.  has  values  along  both  the  positive  and  negative  halves  of  the  imaginary 
axis.  This  causes  a  significant  difference  between  the  leading  edge  and  trailing 
edge  o.l.t.f.  when  the  out  of  phase  aerodynamics  are  included  in  the  binary 
flutter  mechanism  (i.e.  [B]  *  [o]  ).  This  is  shown  on  the  right  hand  side  of 

figure  6. 


It  shows  that  if  the  out  of  phase  control  surface  aerodynamics  are  ignored, 
active  flutter  control  can  be  achieved  using  the  leading  edge  control  driven 
by  displacement  (or  acceleration)  output  at  O/A,  but  only  with  poor  phase  margins. 
On  the  other  hand,  using  the  trailing  edge  control  with  velocity  feedback  from 
O/A  will  give  excellent  phase  margins. 

However,  when  the  control  surface  out  of  phase  aerodynamics  are  included, 
the  leading  edge  control  o.l.t.f.  is  almost  unchanged  but  the  trailing  edge 
control  o.l.t.f.  loses  almost  all  of  its  potential  phase  margin  (ref.  3). 

CONTROL  LAW  DESIGN 


Ref.  6  and  7  show  by  analytical  studies  that  the  stability  magins  specified 
in  section  3  can  often  be  achieved  with  a  simple  control  law  using  a  single 
transducer  and  a  single  control  surface.  This  is  done  by  careful  choice  of  the 
transducer  position.  However,  section  4  reports  that  the  phase  margin  requirements 
(+60°)  of  the  YF-17  Design  Objectives  cannot  be  achieved  using  a  single  transducer 
with  either  the  leading  edge  or  trailing  edge  control  surface  (see  also  ref.  8). 
Therefore  a  multiple  output  control  law  is  required. 

1  Optimal  Control  Theory 

An  analytical  control  law  which  satisfies  the  Design  Objectives,  and  does 
so  using  minimum  control  surface  motion  in  tunnel  turbulence,  can  be  derived  from 
Optimal  Control  Theory.  This  is  fully  described  in  references  6  and  9  and  out¬ 
lined  in  figure  7. 

The  theoretical  derivation  follows  recognition  that  the  Matrix  Riccali 
Equation  for  solving  the  Linear  Optimal  Regulator  Problem  can  be  put  into  Pole 
Placement  form.  Then  letting  the  state  weighting  matrix  fQj  in  the  cost 
function  be  zero  leads  to  the  simple  criterion:  the  control  law  allows  the  closed 
loop  poles  to  be  the  same  as  the  open  loop  poles  except  that  the  real  part  of  the 
unstable  open  loop  poles  change  sign.  This  means  that  the  Nyquist  Plot  will  consist 
only  of  an  anticlockwise  circle  with  its  centre  at  -1,0  and  of  unit  radius,  for 
every  unstable  root.  In  other  words,  the  feedback  signal  contains  only  flutter 
mode  velocity  (fig.  8). 

However,  this  control  law  formally  requires  the  sensing  of  all  the  system 
states  (e.g.  20  states  for  10  second  order  degrees  of  freedom)  and  the  feedback 
signal  to  be  a  weighted  sum  of  these  states.  This  is  obviously  not  practical. 

Fortunately,  for  this  configuration  of  the  YK-17  model  the  flutter  mechanism 
is  a  pure  binary  fairly  well  separated  in  frequency  from  adjacent  roots.  Therefore, 
the  control  law  can  be  designed  for  the  flutter  binary  and  be  expressed  as  the 
weighted  sum  of  the  displacement  and  velocity  output  at  two  transducer  positions 
with  a  correcting  filter  to  allow  for  the>  actuator  transfer  function  (fig.  9). 


Such  a  control  law  was  designed  by  Pole  Placement  computation  for  the  YK-17 
model  at  1.7QF  using  the  output  of  the  two  outboard  transducers.  The  resulting 
Nyquist  Plot  for  the  10  d.of  f.  system  with  this  control  law  is  shown  in  figure  ID. 
It  gives  the  required  circular  anticlockwise  lobe  and  good  stability  margins  at 
all  frequencies. 


A  Practical  Control  Law  Design  Procedure 


However,  a  more  practical  approach  can  be  made  by  direct  use  of  the  known 
desirable  anticlockwise  circular  Nyquist  Plot.  That  is  to  take  the  fotu  component 
o.l.t.f.  (i.e.  velocity  and  displacement  outputs  foi  two  transducer  positions) 
and  find  the  four  gains  required  to  combine  them  to  give  a  least  squat  os  lit  to 
the  perfect  circle.  In  this  way  the  actuator  transfer  function  can  be  included 
in  the  basic  o.l.t.f.,  significant  response  in  adjacent  modes  over  the  flutter 
root  bandwidth  can  be  allowed  for,  and  empirical  o.l.t.f.  can  he  used  if  necessaty. 


Fig.  11  shows  the  results  of  doing  this  using  the  two  outtxiaid  at  t  e  1  ot  ome  t  ei  s 
an<l  including  the  actuator  transfer  function.  Partial  integration  of  the  acielot- 
ometer  outputs  is  used  to  avoid  non  zero  response  at  zero  frequency.  Tin  lesults 
are  good.  The  flutter  lobe  fits  the  desired  circle  well  enough  ami  the  only  sig¬ 
nificant  open  loop  stable  mode  which  is  threatened  by  rinsed  loop  instability  is  Du 
calculated  18  Hz  wing  tip  launcher  rail  fundamental  bonding  modi’.  Howevig  ,  tin 
aerodynamic  damping  in  this  mode  is  only  ().  l'V  critiial  and  stnii  tui.il  damping  (not 
included  in  these  tabulations)  will  enoimously  attenuate  this  loin. 
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A  notch  filter  can  then  provide  any  further  attenuation  necessary. 

Fig.  12  shows  the  stability  margins  between  QF  and  2.0QF  for  the  control 
law  designed  at  1 . 7QF  shown  in  fig.  11.  Those  stability  margins  are  excel!* »t 
and  close  to  those  of  the  Design  Objectives.  The  block  diagram  for  this  control 
law  is  the  upper  horizontal  path  of  fig.  13.  This  will  be  referred  to  as  the 
four  gains,  two  transducers  (4G/2T)  control  law. 

Effect  of  4G/2T  Control  Law  on  Flutter  Aerodynamic  Coefficients 

Any  active  flutter  control  law  can  be  thought  of  as  effectively  causing  a 
change  to  tne  flutter  generalised  force  coeff icient s.  To  investigate  this,  the* 

4G/2T  control  law  was  designed  for  the  flutter  binary  of  the  YF-17  model  at  1.7QF 
using  the  leading  and  trailing  edge  control  surfaces  in  turn. 

Although  all  eight  generalised  aerodynamic  force  coefficients  wen*  effectively 
changed,  the  major  effects  were: 


(a)  leading  edge  control  surfaces 

-  0.2  x  work  done  on  torsion  mode  by  in  phase 


bending  modi'  aerodynamic*: 


out  of  phase 


0.7  x  " 


n  in  phase 


x  u  ii 


out  of  phase 


(b)  trailing  edge  control  surfaces 


-0.3  x  work  done  on  bending  mode  by  in  phase' 


torsion  modi*  aei  odvnami cs 


out  of  phase 


"  in  phase 


bending  mode  aerodynamics 


"  out  of  phast 


This  shows  that  the  leading  and  trail  inn  .  dge  control  surfaces  do  work 
primarily  on  the  torsion  and  bending  modes  i espect i vel y.  The  main  stabilising 
effect  comes  from  changing  the  sign  of  a  cross  aerodynamic  stiffness  while  the 
other  changes  provide  the  required  stability  margins  (ref.  "). 

b.  FURTHER  CONTROL  LAW  DESIGN 

o . 1  Analytical  Studios  using  Later  Data 

The  work  described  so  far  (ref.  v>)  has  been  done  using  the  earlier  analytical 
data  (10  d  of  f  system)  which  included  seven  calculated  structural  modes.  This 
work  was  repeated  (ref.  10)  using  later  data  (l  l  d.  of  I  system)  which  included  in 
GVT  mode  s . 

Fig.  13  shows  the  deflections  over  the  wing  in  these  10  structural  modes. 

The  first  two  are  the  fundamental  bending  and  torsion  modes  of  the  flutter  binary. 

Fig.  14  compares  the  o.l.t.f.  for  input  to  the*  leading  edge  control  sui face 
actuator  and  velocity  output  from  the  outboard  accelerometer  positions  at  1.7QF. 

For  both  sets  of  data  the  potential  phast'  margins  for  single  output  control  law 
are  inadequate. 

b • 2  Extra  Feedback  Filters  Defined  by  Northrop 

From  their  experience  in  earlier  tests  (ref.  1)  Northrop  stated  that  extra 
filters  were  required  to  avoid  closed  loop  instabilities  at  frequencies  well 
separated  from  the  flutter  root  bandwidth.  These  art*  two  notch  filters  (  14  and  ‘>d 
Hz)  and  a  second  order  low  pass  filter  (42  Hz  cutoff)  to  protect  high  frequency 
roots,  anti  two  high  pass  filters  (  r».  3  anti  1 .  <>  Hz  cutoff)  to  avoid  rigid  Uidv  instab¬ 
ility  (fig.  13  and  lb). 

To  study  the  effect  of  these  filters,  as  wel 1  as  the  leading  edge  control 
surface  actuator,  the  single  output  o.l.t.f.  for  the  outlvoaid  accelerometer  positions 
an<l  the  Nvquist  Plots  for  the  4G/2T  control  laws  were  produced  for  the  six  filter 
combi na t ions  tabulated  in  fig.  lb. 


Fig.  1 7  shows  the  attentuation  and  phase  lag  produced  by  each  filter  combin¬ 
ation  at  the  1 . 7QF  flutter  root  frequency. 
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Fig.  18  shows  the  effect  of  these  filter  combinations  on  the  anticlockwise 
lobe  of  the  o.l.t.f.  for  leading  edge  control  surface  input  and  displacement 
output  at  O/A.  Whereas  the  mean  phase  of  the  lobe  varies  with  filter  combination 
as  would  be  expected  from  fig.  17,  the  potential  gain  and  phase  margins  for  a 
single  output  control  law  reduces  dramatically  as  the  number  of  filters  increase** 
such  that  for  FCb  (the  filter  combination  used  eventually  in  the  wind  tunnel  test) 
they  are  only  *1.3  and  +4  . 

Fig.  19  shows  the  effect  of  each  filler  combination  on  the  Mutter  lobe  of 
the  4G/2T  control  law  Nyquist  Plot  at  the  Design  Condition  of  1.7QF.  Shown 
dashed  on  each  flutter  lobe  is  the  ideal  circle  of  unit  radius  and  with  centre 
at  -1,0. 


Although  the  circle  fit  is  not  very  good  after  FC1  and  FC2,the  stability 
margins  are  excellent  for  all  filter  combinations. 

The  four  feedback  gains  (J1  to  J4)  were  calculated  by  the  least  squares 
fit  procedure  described  in  section  5.2  for  each  of  the  six  filter  combinations 
as  shown  in  figure  20. 

Fig.  21  shows  the  closed  loop  flutter  curves  for  the  control  law  f or  FC* 
designed  at  1.7QF.  For  clarity,  only  the  roots  of  the  first  two  structural  modes 
(i.e.  the  flutter  binary)  are  shown.  They  show  no  tendency  towards  closed  loop 
instability.  Also,  a  Nyquist  Plot  study  showed  the  same  for  the  other  roots. 

Section  5,1  suggested  that  the  frequency  of  the  first  two  roots  would  not 
change  as  the*  feedback  gain  increased.  Fig.  22  shows  that  this  is  true  for  the 
flutter  root  but  not  for  the  first  root.  This  is  because  the  filters  prevent  the 
anticlockwise  lobe  for  FCr>  in  fig.  19  from  being  a  perfect  circle  (as  it  is  for 
FC1)  so  that  there  is  significant  feedback  of  first  mode  response. 

Because  it  is  only  the  filter  properties  over  the  flutter  root  frequency 
bandwidth  which  affect  the  value  of  the  four  gains  and  the  flutter  root  stability 
margins,  it  would  be  much  better  to  design  a  bandpass  filter  which  introduces  the 
minimum  of  phase  and  amplitude  change  over  this  bandwidth  while  providing  the 
required  amount  of  attenuation  and  phase  change  at  the  lower  and  higher  frequencies 
where  closed  loop  instabilities  may  be  introduced.  This,  of  course,  would  require 
an  accurate  definition  of  the  attenuation  and  phase  requirements  to  avoid  closed 
loop  instabilities. 

7.  WIND  TUNNEL  TEST  RESULTS 


During  the  wind  tunnel  tests  in  October  1979,  the  AFFDL  Design  Objective  was 
reduced  to  increasing  the  flutter  dynamic  pressure  by  70%  without  any  requirement 
to  demonstrate  stability  margins. 

It  had  been  expected  that  the  analytical  control  law  would  not  he  adequate 
to  provide  the  required  stability  margins  in  the  test.  It  had  therefore  been 
planned  to  obtain  empirical  o.l.t.f.  and  use  the  least  squares  fit  procedure*  on 
th«  se  to  produce  empirical  4C./2T  control  laws. 

However,  the  70%  increase  in  flutter  dynamic  pressure  was  achieved  using  the 
law  calculated  using  the  13  d.  of  f.  analytical  data  for  the  additional  filters 
that  Northrop  considered  necessary  i.e.  filter  combination  b  (fig.  lb). 

As  the  tests  proceeded  in  dynamic'  pressure,  ihc  stability  margins  were 
partially  checked  by  reducing  the  control  1 ^w  •  »v«  rail  gain  down  to  0.7  of  its 
nominal  value  and  also  introducing  up  to  30  phase  lag  at  the  flutter  root  frequency 
by  means  of  an  all  pass  filter.  Resulting  from  this,  the  control  law  was  tested 
at  1 . 7QF  with  only  0,9  and  0,8  of  the  nominal  gain  and  with  no  phase  lag.  This 
was  because  the  tests  at  lower  dynamic  pressure*  had  predicted  that  at  1 . ?0F s  the 
flutter  root  would  not  be  stable?  for  0.7  of  nominal  gain;  and  for  nominal  gain, 
with  or  without  any  phase  lag,  the  model  response  to  tunnel  turbulence  in  the  IV Ion 
Side  Bending  Mode  at  about  13  Hz  would  cause  motion  of  the  leading  edge  control 
surface  which  could  damage  the  actuator. 

Because  of  this  response  to  tunnel  turbulence,  it  was  n«»t  possible  t"  obtain 
good  empirical  o.l.t.f.  alxive  1 . 2QF  (no  p.s.f.).  Howrvi  i  ,  from  the  comparison  of 
the  analytical  and  empirical  o.l.t.f.  at  1 . 2QF  it  was  possible  to  obtain  empirical 
correction  factors  for  the  analytical  data  concerning  the  flutter  and  1  *•  Hz  lobes 
(ref.  11). 

The  effect  of  applying  these  correction  factors  t«>  tin’  analytical  Nyquist  Plot 
at  1 . 7yh  is  shown  in  fig.  21,  The  flutter  lobe  is  reduced  by  a  factor  of  0.7;  the 
1  t  Hz  lobe  is  increased  by  a  factor  of  2  and  has  a  b()  phase  lead  relative  to  the 
negative  real  axis. 

This  c  harnir  to  the  analytical  Nyquist  Plot  makes  it  consistent  with  the  test 
results  which  limited  the  1 , 'QF  tests  to  0,9  and  0.8  times  the-  nominal  gain. 


Given  this  information,  it  appears  that  good  empirical  o.l.t.l.  could  be 
obtained  at  1 . 7QF  by  merely  introducing  a  notch  filter  into  the  contiol  Jaw  to 
prevent  the  large  control  surface  response  to  tunnel  turbulence.  Then  those 
o.l.t.f.  could  be  used  to  design  a  control  law  with  good  stability  margins. 

CONCLUDING  REMARKS 

The  YF-17  model  tested  had  extreme  flutter  characteristics  which  caused  a 
new  control  law  design  procedure  to  be  developed  to  obtain  good  staid  lity  margins. 
This  procedure  gave  a  control  law,  based  upon  analytical  data  and  a  leading  edge 
control  surface,  which  was  robust  enough  to  achieve  the  required  7 03?  increase  in 
flutter  dynamic  pressure. 

It  could  be  that  when  designing  for  active  flutter  control  of  an  aircraft, 
the  analytical  data  is  not  adequate  so  that  good  quality  empirical  data  may  be 
necessary.  Also,  only  a  trailing  edge  control  surface  may  be  available. 

Therefore,  the  following  work  on  the  YF-17  model  would  be  useful, 

(1)  Analytical  studies  to  find  the  relative  effectiveness  of  U*  leading  and 
trailing  edge  control  surfaces  for  active*  flutter  control. 

(2)  Wind  tunnel  tests  to  achieve  1,7QF  using  a  control  law  based  upon  analytical 
data  and  the  trailing  edge  control  surface. 

(3)  Wind  tunnel  tests  designed  to  obtain  good  quality  open  loop  transfer  functions 
by  separately  exciting  with  the  leading  and  trailing  edge  control  surfaces 

up  to  1.7QF. 

(4)  Wind  tunnel  tests  of  control  laws  based  upon  the  empirical  data  of  (  i)  above, 
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